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1. 

AN  HDCP^IiraTTAL  INVESTIGATION  0?   THE  G-AINS   OF   DIELECTRIC 

ANTENNAS 

INTRODUCTION 

The  idea  of  using  dielectric  materials  for  the 
control  of  the  direction  of  electroma.^netic  i^adiation 
V7as  probably  first  conceived  by  Hertz  in  the  Nine- 
teenth Century  ,  when  he  observed  the  strong  influence 
of  pitch  blocl^s  on  this  little  ^onderstood  form  of 
energy  flov;.   Subseauently,  in  the  early  part  of  the 
Ti-jentieth  Century  a  fev  investigations  were  reported 
concerning  the  use  of  dielectrics  as  vrave  guides, 
but  there  vere  only  occasional  contributions  to  this 
theory.   Rapid  progress  in  this  direction  v.'as  not 
begun  until  the  availability  of  microwave  equipment 
in  the  1930 *s.  By  1939,  dielectric  vave  guide  theory 
had  been  thoroughly  studied  and  a  few  experiments 
had  been  attempted  in  the  use  of  dielectrics  as 
directional  ra,diators.  'vith  the  stimulus  of  the  v;ar, 
and  the  ranidly  increasing  need  for  more  highly 
directive  antennas,  dielectric  rods  were  studied, 

1.  Superscript  nujnbers  refer  to  references  listen  in 
Bibliography, 
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tried,  and  put  Into  use  beglnnins  in  the  early  1940' s. 
In  1943  the  U.S.  Navy  put  into  service  a  radar 
equir)ment,  designed  by  the  Bell  Telephone  Laboratories, 
V7hich  employed  an  array  of  forty-two  rods  of 
Polystyrene  to  give  a  sharp  directional  beam  in  the 
horizontal  plane.   (Radar  Equipment  Mark  Till). 

The  G-ermans,  about  the  same  time,  began  to  use  a 

2 
directional  antenna  on  aircraft  which  consisted  of 

four  dielectric  rods  arranged  to  sive  a  symmetrical 

pattern.   Subsequently,  a  large  number  of  dielectric 

antennas  have  been  built,  and  an  ever  increasing 

demand  for  these  small  size,  highly  directional 

devices  has  encouraged  more  and  more  a-ttention  to 

their  theory  of  operation. 

It  should  be  pointed  out  in  the  beginning  that 

y     the  classical  methods  of  antenna  design,  based  on 

Mar\/ell's  equations,  ha,ve  not  been  successfully 

applied  to  dielectric  antennas  whose  dimensions  liave 

been  of  the  order  of  a  wavelength.   As  clearly  Dointed 

out  by  Stratton-^  in  his  book,  "Slectromagnetic  Theory", 

Maxwell's  equations  are  valid,  in  general,  only  to 

large  scale  phenomena,  and  attempts  to  extend  them 

to  the  so  called  "Quasi-microscopic  state"  lead  only 

to  partial  success.  Although  many  of  the  laws  of 
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clasr.ical  electrodyTiamlcs  apply  to  submicroscopic 
domains,  It  is  not  possible  to  apply  them  seneral'ly. 
On  the  other  side,  thei-'e  is.  quantum  electrodynarrlcs, 
v;hich  at  this  x-zritin.'^,  have  not  been  extended 
throuf?h  th3  transition  areas.   The  statistical 
avera-^-es,  over  a  lar.^e  number  of  atoms,  of  the 
equations  of  quantum  electrodynamics  ultimately 
must  lead  to  the  relations  derived  by  Kaxwell. 

Design  of  dielectric  antennas  to  date  has  been 
based  largely  on  cut  and  try  methods.   Some  of  their 
characteristics  have  been  determined  by  various 
experimenters,  and  considerable  theoretical  vv-ork 
has  been  done  in  an  endeavor  to  exp'lain  their 
operation.  But  there  is  no  mass  of  data  for  the 
theorist  to  draw  upon,  nor  is  there  an  adeouate  theory 
at  the  present  time  to  enable  calculation  or  prediction 
of  the  characteristics  of  a  given  dielectric  shape. 
It  is  the  purpose  of  the  -rriter,  in  undertaking,  this 
work,  to  determine  experimentally  the  "concentrating 
action"  of  certain  of  these  antennas,  and  to  express 
this  in  terms  of  nain  relative  to  an  isotropic 
oscillator;  or,  './hat  amounts  to  the  sa-me  thinr;,  to 
express  this  gain  as  a  ratio  of  effective  area  to 
actual  area.   It  is  believed  that  these  data,  thus 
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obtained,  may  be  of  use  In  subsequent  invest! r^atlons, 
or  for  other  investigators,  to  ultimately  develop 
a  satisfactory  theory  of  the  function  of  the  dielectric 
antenna.   Finally,  available  theoretical  approaches 
will  be  discussed  objectively  in  the  lir^ht  of  the 
experimental  results  obtained  herein. 

'•■'ork  of  other  investi.^atorsi 

Mallach  ,  in  Q-ermany,  presented  a  thesis  in 
1939  wliich  reported  considerable  work  in  the  use  of 
dielectrics  to  produce  a  directional  beam.  He  also 
produced  a  simple  theory  v^hich  explains  oualiltatiT^^Jy 
the  action  of  the  dielectric,  but  failed  quantitatively 
to  verify  all  of  the  knovm  data.   Taking  a  dielectric 
filled  vrave  guide,  whose  walls  are  considered 
perfectly  conducting,  excited  at  one  end  and  open  at 
the  other,  it  is  well  kno^^ni  that  radiation  is  emitted 
from  the  open  end,  and  this  radiation  is  directional 
alonp;  the  projected  axis  of  the  guide.   If  the 
conducting  sheath  is  removed  from  a  part  of  the  guide, 
nearest  the  open  end,  leaving  some  of  the  dielectric 
rod  protruding,  radiation  will  occur  both  from  the 
end  of  the  dielectric  and  from  its  side  v-alls  as 
well.  Mallach* s  theory  is  based  on  the  consideration 
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that  each  olement  of  the  surface  acts  as  a  radiator 
and  the  phar.e  of  this  excitation  depends  on  the 
extended  len-^th  of  the  rod..  He  assumes  in  this  work 
that  the  v;ave  is  linearly  solarized  and  of  constant 
intensity  and  Is  transmitted  dovm  the  rod  without 
attenuation  or  reflection.   Furthermore,  he  does  not 
consider  the  rod  cross  section  area  except  for  the 
determination  of  the  wavelength.   He  also  considers 
only  small  angles  of  diffraction,  takin'^  no  account 
of  the  spherical  nature  of  the  wave.   Mallach 
reported  a  considerable  amount  of  experimental  data 
on  both  rods  and  dielectric  tubes,  and  arrived  at 
some  sim;ole  formulas  for  choosing  dimensions.   For 
the  sizes  of  rods  and  tubes  that  he  used,  his  simple 
theory  and  formulas  are  satisfactory,  but  are  of  no 
use  in  obtalnln-s  lobe  patterns  for  the  entire 
radiation  field  nor  for  computing  gains. 

•  Southworth' ,  at  the  Bell  Telephone  Laboratories, 
vrhose  early  publications  on  wave  guides  and  v^ave 
guide  theory  are  v.'ell  known,  has  done  considerable 
experimenting  for  more  than  a  decade  i;ith  tapered 
dielectric  rods,  both  rectangular  and  circular,  of 
diameters  of  a  half  wavelength  or  less.  He  and  his 
associates  have  experimentally  established  that 
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the  P^aln   decreases  as  the  attenuation  in  the  rod 
increases. 

o 

v/ilkes  ,  in  a  thesis  presented  to  the  University 
of  Paris  in  1946,  examined  the  problem  from  the 
point  of  vieK  of  the  lens  action  of  the  dielectric 
on  radiation  of  wavelengths  comparable  to  the 
dimensions  of  the  lens.   He  pointed  out  that  the 
apparent  index  of  refraction  plays  an  important 
role  in  the  energy  gain  of  the  dielectric  radiator. 
By  Dlotting  experimentally  obtained  energy  gains 
versus  length  of  exposed  dielectric  he  has  shovrn  a 
correspondence  of  periodic  gain  to  Snell^s  law  of 
refraction.  Furthermore  he  has  shown  that  tho 
apparent  index  of  refraction  is  itself  not  constant, 
but  depends  on  the  exposed  length  of  dielectric  and 
a  cross  section  dimension.   A  formula  for  calculating 
the  apparent  index  from  the  real  index  of  refraction 
is  developed  in  his  paper.  V/'ilkes  also  computed  some 
radiation  patterns  based  on  analogy  to  more 
conventional  antennas,  v/hich  verified  very  nicely  the 
experimental  data  for  the  major  lobe.   But,  as  his 
paper  points  out,  his  assumptions  used  would  cause 
considerable  errors  in  the  side  lobes.   This 
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calculation  assumes  constant  wavelength,  vrhich  he 
sho'w's  to  vary  somewhat  throughout  the  lens  length. 
He  also  assumes,  for  tills  calculation,  constant 
displacement  currents.   The  complexity  of  an  exact 
formulation  of  wavelength  and  intensity  would  lead 
to  expressions  v;hich  v/ould  be  extremely  difficult 

if  not  impossible  to  solve. 

11    1  o    1  f^ 
Several  investigators  *   '    have  subsequently 

published  attempts  at  calculation  of  radiation 

patterns.  All  of  the  recent  attemrtts  knovm  to  the 

author  are  based  on  the  formal  solution  of  the  field 

equations  given  in  Schelkunof f  *  s  well  knov.Ti  book^, 

as  applied  to  graves  In  dielectric  V7ires,  and  expressed 

In  terms  of  cylinder  functions.   Elsasser^  in  1949 

published  an  analysis  of  the  various  non- radiative 

modes  in  a  circular  dielectric  rod,  in  vrhich  he  also 

pointed  out  a  perturbation  method  of  calculating 

dielectric  loss,   Blsasser,  In  this  paper,  did  not 

attempt  calculation  of  radiation  from  a  finite  length 

rod.   This  theoretical  study  is  in  substantial  ■ 

agreement  v/lth  an  experimental  investigation  of 

dielectric  cylinders  as  v;ave  guides  performed  by 

Chandler   at  RCA  laboratories. 
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Marston   of  the  U.S.  Naval  Research  Laboratory 

has  attempted  calculation  of  radiation  from  dielectric 

rods  of  Tiniform  circular  cross  section.  He  employed 

the  method  of  assumed  fields  based  on  Schelkunoff * s 

1  p 
Equivalence  theorems  ,  published  in  the  Bell  System 

Technical  Journal  in  1936.  Marston  computed  the 

radiation  field  arisin^^  from  equivalent  sources  at 

the  sides  of  the  rod,  and  separately  for  those  at 

the  ends  of  the  rod,  and  points  out  that  maximum 

gain  in  the  axial  direction  is  achieved  v:hen  the 

total  radiation  from  the  end  is  in  phase  v/ith  the 

radiation  from  sides  (for  modes  of  type  n  =  1  only). 

Marston' s  calculations  also  show  that  the  actual 

values  of  gain  for  dielectric  circular  rods  cannot 

in  general  be  predicted  from,  ordinary  array  theory, 

but  prospects  of  achieving  this  exist  through  the  use 

of  his  complex  expressions  for  the  components  of 

electric  field  intensity.   However,  all  of  Marston 's 

assumptions  are  not  in  agreement  with  findings  of 

other  investigators;  for  example  he  assujnes  constant 

vravelength  within  the  dielectric  body.   If  this 

variable  v;ere  substituted  for  the  constant,  it  would 

make  his  equations  almost  hopelessly  complex. 
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A  great  deal  of  the  knovm  Information  on 

1? 
dielectric  antennas  has  been  assembled  by  I-lcKinney 

in  a  dissertation  presented  to  the  University  of 

Texas  in  June,  1950.  He  re-norted  also  a  very 

substantial  a-raount  of  experimental  data,  as  well  as 

a  3lir;htly  different  method  of  attempted  calculation 

of  the  fields.   The  method  presented  is  also  ba-sed 

on  Schelkunoff 's  Equivalence  Theorems, .  but  uses  an 

extension  devisad  by  Horton,  in  v.'hich  he  computes 

fictitious  electric  and  magnetic  current  sheets 

about  the  dielectric  vhich  are  a  useful  aid  in  the 

analysis.   Then  he  obtains  the  vector  potentials  at 

a  point  In  space  due  to  these  sheets  and  computes  the 

E  and  H  fields  at  the  point.   In  his  calculations  he 

neglects  the  radiation  from  the  end  of  the  rod, 

neglects  reflections  In  the  rod,  and  neglects  any 

changes  which  mi^ht  occur  to  the  v.'avelen^th  of  the 

wave  v.'ithin  the  dielectric.   McKlnney's  experimental 

data  agreed  v/lth  his  theoretical  calculations  quite 

closely.   For  the  sizes  and  shapes  used  and  for  the 

TMq-,  mode,  it  appears  that  his  solutions  are  quite 

precise.  Hov;ever,  his  report  does  not  show  any 

calculations  of  patterns  for  dielectric  antennas 

excited  in  the  domins^nt  mode,  which  is  the  mode  most 
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commonly  used,  nor  does  his  theory  conta-in  a  precise 
statement  of  the  optimum  j^alns  v/hich  ras-y  be  obtained. 
The  TMq,  mode  gives  a  null  on  axis,  which  Is 
unacceptable  for  most  cases  where  dielectric  antennas 
are  required. 

It  should  be  pointed  out  that  no  theoretical 
treatment  of  any  dielectric  antennas  other  than  those 
of  circular  cross  section  have  been  found  in  any  of 
the  literature.  As  v:ill  be  brought  out  later  in  this 
essay,  an  attempt  at  a  straightf oirvjard  solution  to 
the  rectangular  shaped  antenna  by  any  of  the  classical 
methods  would  lead  to  exceedin?3ly  complex  eouations 
v^hose  solutions  v:ould  involve  serious  difficulties. 
Nevertheless,  in  a  great  many  applications,  the 
rectangular  antenna  is  the  one  preferred  due  to  the 
predominant  use  of  rectangular  wave  guide  which 
assures  no  rotation  of  the  plane  of  polarization. 
Hence  this  essay  includes  experimental  gains  obtained 
for  both  rectangular  and  circular  cross  section 
guides,  including  such  correlation  to  existing 
theories  as  are  possible. 
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CHAPTER  1. 

RliJVlSl-/  0?  BASIC  LAV'S  OF  ELECTROFAdNETIC  RADIATION 

The  objective  of  Chapter  1  Is  to  develop  the 
eouations  for  the  propagation  of  electromap;netic 
vraves  In  a  circular  dielectric  rod,  infinitely 
long,  imbedded  in  a  homofjeneous  medium.   The  starting 
point  will  be  Max^jell's  equations.   The  development 
will  be  carried  out  for  the  cylindrical  case  because 
treatments  of  this  are  available  in  the  literature. 
Throughout  this  work  the  I-IKS  system  of  units  will 
be  used. 

It  has  been  pointed  out  in  the  introduction 
that  MaxT-:ell*s  equations  hold  with  certainty  only 
for  macroscopic  phenomena.   Accordingly,  we  make 
the  assumption  throughout  Chapters  1  and  2  that  only 
macroscopic  conditions  are  considered,  that  the 
media  considered  are  isotropic  and  homogeneous,    > 
and  that,  except  for  stated  discontinuities,  the 
physical  properties  in  the  neighborhood  of  every 
point  are  continuous.   Under  these  conditions,  the 
following  Maxv;ell*s  equations  have  been  firmly 
established. 


(1)   yx£  t||=  O 
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<2'       vx/y-^  =  ^ 

where  the  vector  quantities  E,  B,  H,  D  and  J  are 
respectively  the  electric  intensity,  the  magnetic 
induction,  the  magnetic  intensity,  the  electric 
displacement  and  the  current  density. 

Since  the  divergence  of  the  curl  of  any  vector 
Is  identically  zero,  we  obtain  by  taking  the 
divergence  of  (1) 

(3)  v-|f-=o 

and  under  the  condition  of  continuity  imposed  above 

Since  the  divergence  of  the  magnetic  intensity  does 
not  change  with  time,  and  certainly  at  some  time 
it  did  not  exist,  it  follows  that 

Likewise,  taking  the  divergence  of  (2)  and 
commuting  the  partial  and  del  signs 


'^'   V-  tT^  ^v-  d-^o 
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Nov^  If  vre  denote  the  charge  density  at  any 
point  by  f ,  the  eouatlon  of  continuity  or  the 
conservation  of  charge  states  that  the  time  rate  of 
change  of  f plus  the  divergence  of  the  current  density 
is  identically  zero,  or 

(7)  V  •  J  ^  H  =  (9 

Substitutins  for  V-j"  in  (6) 


or 


(9)   |7(V-D-fj=0 


Applyln,;-  a,5ain  the  argument  that  if  the  time  rate 
of  change  of  the  quantity  is  zero,  the  quantity 
itself  must  be  zero  if  it  ever  has  vanished,  it 
follows  that 

(10)   V-  D  -  f 

This  tells  us  that  the  charges-  distributed  vjith  the 
density p  constitute  the  sources  of  the  vector  D, 

Electric  and  Ma,cr,netic  Polarization 

In  order  to  describe  the  electromagnetic  state 
of  a  substance,  it  is  convenient  to  define  tv:o 
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additional  vectors,  the  electric  and  the  magnetic 
polarization  vectors.   By  deflnl  ion 

(11)   P  =D  -  g.£ 

(12)   M  --  i;fi  -H 

where  ,£„  is  the  inductive  ca-nacity  in  free  space 

and  Mo    is  the  permeability  of  free  space.   Therefore 

in  free  space,  the  polarization  vectors  vanish. 

Substituting  the  values  for  D  and  H  from  (ll) 
and  (12)  into  (2)  and  (lO),  v/e  obtain  the  follov/ing 

(13)  V^^B-^oio   II  ^^o(jf|^  tVx/Vl) 


(14)   V  E  '-  iC?'V-  P) 


In  any  physical  media  \fhlch  is  isotronlc  and 
In  which  the  field  vectors  are  linear,  the  Electric 
DisiDlacement  and  the  Magnetic  Induction  are  related 
to  the  corresponding  field  intensity  by  the  constants 
of  the  medium 

(15)  D  ^  £E 
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(16)  B  -  ^  H 

From  the  above  enuations  (13-lA)  it  Is  possible 

to  accoiant  for  material  bodies  in  an  electromagnetic 

field  by -an  equivalent  charge  of  density -VP   and 

an  equivalent  distribution  of  current  density  ^  ■'■  V-^/M 

It  is  convenient  at  this  point  to  define  two 
additional  coefficients  of  proportionality.   The 
electric  and  magnetic  susceptibilities  of  any  media, 
excluding  ferromagnetic  material,  will  be  denoted 
by^Q  ^'^^^m  ^^'^   ^^®  defined  by 

(17)  P--X,LE 

(18)  M=X«,H 

We  can  define  these  scalar  quantities  in  this  vray, 
since  in  isotropic  media,  the  polarization  vectors 
are  parallel  and  proportional  to  the  corresponding 
fiteld  vectors. 

It  is  also  convenient  to  define  the  relative 
dielectric  constant 


(19)      e  -'  f 

and  the  relative  permeability 


(20)     K.   --   I; 
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Prom  (11),    (15),    (17).    and  (19)   it  follows  that 

(21)    X,  -=  e -I 

and  from  (12),  (16),  (18),  and  (20)  that 

(22)  'X^    ^    K^-l 

Vector  and  Scalar  Potentials 

To  facilitate  the  analysis  of  electromagnetic 
fields  it  is  useful  to  define  some  auxiliary  functions 
called  potentials.   These  may  be  shown  as  follows: 

Equation  (5)  shows  that  the  divergence  of  the 
Magnetic  Induction  Vector,  B,  is  zero  at  all  times. 
According  to  this,  the  field  of  B  is  al^^ays  solenoidal, 
and  therefore  B  can  be  represented  as  the  curl  of 
another  vector  A, 

(23)  B  =  yXA 

Substituting  this  value  of  B  into  equation  (l) 
gives 

3  A\    _ 


(24)  V^  {E^    U)    ==  0 


Thus  the  field  of  the  vector  (E  f  |j  )  is 
irrotational  and  equal  to  the  gradient  of  a  scalar 
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function  '.fhlch.  v;e  may  call  (j), 

(25)  Sf-^  =  -V^      or  E  =  -V^-|^ 

The  function  A  is  called  the  vector  potential  of 
the  field,  and  the  (^  Is  the  scalar  potential. 

It  should  be  noted,  however,  that  the  vector  A 
and  the  scalar  (j)  are  not  unique.  Thus  B  could  have 
been  represented  by  the  curl  of  another  vector 

(26)  B  ==    V^Ao 
where 

(27)  k=Ao-Vf 

where  j//  Is  any  arbitrary  scalar  function  of  position, 
It  would  similarly  follow  that 

(28)  S  ==  -7^^  -  1^ 

and  ^  and  ^  would  be  related  by 

(29)    *  =  i  f  i^ 

Likev/ise  an  infinite  number  of  potentials  describing 
the  same  field  could  be  derived  in  the  same  manner. 

By  the  definition  of  ^,  equation:?  (15),  and 
using  the  above  defined  potential,  equation  (25), 
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(30)  D  =  -S(V(j>   +   If) 

And  by  eouatlons  (l6)  and  (23) 

(31)  H  ^    ivx/\ 
Substltutln-^  (30)  and  (3I)  into  eouation  (2) 


and  multiplying  the  above  hj  M   and  commuting -|r: 
v/ith  V  ^ives 


(32)  ^^^^^MSV^^Mi^^.    Mxl 

By  substituting  (30)  into  equation  (lO)  and  dividing   ' 
by  -  ,5 

(33)  V-V(^  f  V-|^  =-^  P 

^  Now  all  particular  solutions  of  (32)  and  (33) 

lead  to  the  same  electromagnetic  field  v;hen  subjected 
to  the  same  boundary  conditions.   They  differ  among 
themselves  only  by  the  arbitrary  function  (j^   ,      In  order 
to  sT?ecify  the  functions  more  precisely,  the 
condition  may  be  imposed  that  A  and  <^  also  satisfy 
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(34)  V-/\  ^^£^   ^0 

In   order  to  do  this,  It  is  necessary  that  (p   satisfy 

(35)  Vy  -^f  0^  =  V-Ao  ^m£  ^ 


V/here  A  and  (j)     are  particular  solutions  of  (32) 
and  (33). 

The  potentials  A  and  (^  "by  the  aid  of  the  above 
are  completely  specified.   They  are  solutions  of 
the  follov7in!3  equations,  v/hich  are  the  same  as  (32) 
and  (33)  except  for  the  substitution  of  7/l=->v^||- 
from  (34). 

(36)  V^V^A  -  V7'A-^^S^.  =M^ 

(37)  V^   -Mit^     --if 
By  use  of  the  vector  identity 

(38)  V^V>^/\  ^  VV'/)-\7-7/\ 

equation  (36)  reduces  to 

(39)  V^/\-Mt^,-^-M^ 
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The  relations  (23)  and  (25)  for  the  vectors 
B  and  E  are  not  -^^eneral.   To  them  may  be  added  any 
particular  solution  of  the  honogeneous  equations: 


(^l   V^E  f   ^    =  0 


(41) 


Vx  H  -  41  =  0 


(42)   X7B=  0 
(*3)   '7'  D  =  0 

From  the  symmetry  of  this  system  tt  is  evident  that 
it  can  be  satisfied  identically  by 


(44)     D  =  -  V^/\ 


* 


(45)     H  =  -  V^"^- 


from  which  may  be  obtained 


(46)     S  =  -  fV^A"" 


(47)     B  «  'M(v(f>U   if) 
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But  the  new  potentials, A^-*  and  (^^\    must  be  subject  to 
the  conditions 


(48)    y^/l^  -^e0l    -0 


(49)    V'^*  -Mi   0a'  =  0 


(50)  VA'-iME^^O 

A  general  solution  of  the  inhomogeneous  system 
given  by  equations  (l) ,  (2),  (5)  and  (lO)  is,  then 

(51)  B«  VAA  -^  ^  -M  V(p'^ 

(52)  ^   E  «  -  V^  -  1^  -  -^7X4* 

The  functions  A  and  9  are  potentials  of  a  source 
distribution  entirely  external  to  the  region 
considered.   Usually  A  and  (()  are  set  equal  to 
zero,  and  the  potentials  of  all  charges,  vjherever  the 
source,  are  represented  by  A  and  <)). 

At  any  -ooint  where  the  charge  and  current 
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densities  are  zero,  a  possible  combination  of 
potentialsMo  -   0  and  0^  =  O.-  The  function^  is 
then  according  to  (55)  any  solution  of  the  homogeneous 
eauation 


(53)   Vy  -    mS^   =0 


Since  (^Q  =  0  satisfies  this  same  equation  at  this 
point^^raay  be  chosen  in  such  a  way  that  j>   vanishes. 
In  this  case  the  field  can  be  expressed  in  terms  of 
only  the  vector  potential.   Equations  (23)  and  (25) 
become 

(54)  B  =   VM 

(55)  S  «  '  JJf. 
and  (59)  and  (54)  become 


(56)   '7''A-mSjp=0 


(57)    VA   '0 

The  foregoing  paragra-nhs  have  shown  hov;  an  electro- 
magnetic field' in  a  homogeneous  medium  may  be 
reduced  to  a  vector  and  a  scalar  Dotential  v/hich 
will  satisfy  Max^fell's  equations. 
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Hertz  Vectors  or  Polarization  Potentials 

It  is  possible  under  ordinary  conditions  to 
define  an  electromagnetic  field  in  a  homogeneous 
medium  by  means  of  a  vector  fume ti on  alone. 

Considerin.^  at  first,  no  free  cliarges  or 
conduction  currents  within  the  medium,  the  field 
equations  reduce  to  the  homogeneous  form  (equations 
(4o)  through  (43)  ),  Assume  a  vector,  II,  v'hooe 
time  derivative  is  r)ro"oortional  to  the  vecto-r 
potential  A 

(58)  ^  =  ^^  If 

and  substitute  in  (23)  and  (25) 

(59)  B  «  ^£  V/  If 


(60)       E=   -V^-^^l^ 


Using  the  relations  B  -■■u./ia.nd.   D  =  ^F,  substitute  the 
above  values  in  the  homogeneous  equation  (4l).  After 
comjnutlng  the -^  and  V  signs,  the  result  is 

(61)  jj  (VaV^U  tV(/}  i-^S0)  ^  0 
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since  the  scalar  function  (^  is  arbitrary  as  long 
as  it  satisfies  equation  {j>^) ,    it  may  now  be  chosen 
to  also  satisfy 

(62)    (^  =  -  V-E 

If  the  above  chosen  value  of  (^  is  substituted  in 
equation  (61)  and  the  result  integra-ted  v/ith  respect 
to  time,  the  follov/ing  is  obtained 

(63)  VKV)^!  -  vvj    -^MS-^7i   *  a  constant 

The  value  of  the  constant  does  not  affect  the 
determination  of  the  field,  and  so  it  may  be  replaced 
by  zero.   Every  solution  of  the  vector  eouation 


(64) 


V^V^IT  -  VV'E  f^(f  1^  ''0 


determined  an  electric  field  through 


(65)  B  :r^f  VA-ff 


and 


,3 


ax 


{66)  E  =     V  V-I-mS   ^Ji 


Equation  (65)  is  identical  to  (59)  and  {66)    is  the 
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same  as  (60)  ■vrith  the  substitution  of  W'JT        ^^^ 
its  enual  -70  by  virtue  of  (62). 
Applying  the  vector  identity,  equation  (38)>  to 
eouatlon  (64)  gives 

(67)  V  VI  -Mi  0.   =^0 

The  reouirement  that  (j)   satisfy  equation  (53)  is 
met  in  equation  (64)  and  (6?). 

Inasmuch  as  the  Electric  Displacement  Vector, 
D,  as  v.^ell  as  the  Ife-gnetic  Induction  is  solenoidal 
in  a  charge  free  re-^ion,  an  alternative  solution 
could  be  similarly  constructed  of  the  form 


*     >^  aj-*^ 


(68)    A"  -  ^S  f^- 


^  -  -if 


(69)  (/>  =  -  V-I 

(70)  D  ^  ^^£V^§^ 


^nif- 


(71)     H  =  vv-iL^'^e  ^ 

where  II  "  is  any  solution  of  (64)  or  (67), 
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It  therefore  follov:s  that  the  electromagnetic 
field  throughout  a  homogeneous  isotropic  region  in 
v/hich  no  charge  or  current  exist  could  be  resolved 
into  two  partial  fields,  one  derived  from  the 
vector  II  and  the  other  from  II  •  The  source  of 
these  fields  is  exterior  to  the  region.   In  order  to 
shov7  the  physical  significance  of  the  Hertz  vectors 
it  is  necessary  to  relate  them  to  their  sources. 
This  may  be  done  by  finding  Inhomogeneous  ecuations 
from  v;hich  (6?)  is  derived. 

The  vector  D  may  be  expressed  in  terras  of  the 
electric  field  Intensity  and  the  electric  polarization 
vectors,  using  equation  (11),  as  D  ^i^   E  f  P.   Then 
substitute  tb-is  in  the  homogeneous  field  eouatlon. 
Equation  (41)  becomes 


(72)  V'H-<?o|f=|^ 

and  (43)  becomes 

(73)  V-E  '-  -i  v-p 

The  above   tv/o  equations  as  v/ell  as   (4o)    and  (42) 
are   still   Identically   satisfied  hj  (65)    and  {66)  ^ 
provided^  is  re;olaced  by  ^     and  that  II  is  now  any 


27. 


solution  of 


(7A)  7-vr  -m£,^  --ip 


The  source  of  the  vector  II  and  the  electromasnetic 
field  vrhich  it  prescribes  is  a  distribution  of 
electric  polarization  P, 

The  vector  P  is  the  electric  dipole  moment  per  unit 
volume  of  the  medium.   The  vector  II  is  associated 
vrith  a  distribution  of  electric  dipoles  and  may  be 
called  the  electric  polarization  potential. 

Similarly  the  field  associated  ^^.'ith  II  is 
established  by  a  distribution  of  magnetic  polarization. 
The  magnetic  induction  vector,  B,  is  related  to  the 
magnetic  field  intensity,  H,  according  to  eouation? 
(12),  by  B  =^Q  (H  f  M).   Substituting  this  into 
(40)  and  (42)  gives 


(75)   V^E^^.H  --A|^ 
and 

(76)     V-H  =  -  V-M 

If  it  Is  further  stipulated  that  II  shall  be  a 
solution  of 
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(77)   v-vJ*-M,e^  =  -'U\ 

and  VclB.%  yt{    reolaceX  in  eouations  (70)  and  (71), 
o   " 

then  enuatl'ons  (75),  (7^),  (4l)  and  (43)  are 

Identically  satisfied  by  (70)  and  (71). 

The  electric  polarization  P  may  be  induced  in 
the  dielectric  by  the  field  E,  but  it  may  also 
contain  a  part  v/hose  mar^nitude  is  governed  by 
entirely  external  conditions.   Insofar  as  dielectric 
antennas  are  concerned,  only  the  independent  part,  P  , 
Is  of  interest,  v:hich  represents  the  electric  moment 
of  dipole  oscillation  activated  by  external  power 
sources.   This  is  also  true  of  magnetic  polarization. 
To  represent  these  conditions,  equations  (ll)  and  (12), 
expressed  in  modified  form  are 

(78)  E  =  -riO'To) 

(79)  H  =  i  ^  -A^o 

in  which  P  and  M.  are  Dre scribed  and  inde-nendent 
o      o      - 

of  E  and  H  and  where  the  induced  polarizations  of  the 
medium  have  been  absorbed  into  the  parameters  ^   and 
Ui  ,   The  electromagnetic  field  due  to  these  distri- 
butions of  pQ  and  Mq  are  given  by 
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(80)     E  =  VV-H-^g^,    -mV^   I^ 


^X  ,  vTV7.F^_  ../t^. 


(81)    H  =  £  Va|^  ^  VV-J^-^/ff-^ 


Equation  (8o)  above  was  formed,  by  substituting;  in 
eouation  (78)  the  value  of  D  ^iven  by  (70)  and  P 
given  by  (7^),  replacing^Q  and  ^q  by  ^  and  ^. 
Equation  (81)  vas  formed  by  substituting  into  (79) 
the  value  of  B  given  by  (65)  and  M^  frora  (77). 
For  these  eouations  to  be  valid,  however,  it  must 
also  be  stipulated  that  II  and  II  be  solutions  of 
equations  corresponding  to  (74) 


5'X 


(82)   V^J-^(f  5f.  -^Z'. 


(83)   V'  VI* -^^0-!  =-A/„ 


Using  the  vector  identity,  equation  (38),  and 
substituting  (83),  it  iisalso  possible  to  write  (81)  as 
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(84)        H  =  ev  x§f  -t-  V^V^JL^ -Alo 

The  last  relation  shoi-fs  that,  since  B  =  VXA^    t^Xe 
vector  potential  A  can  be  derived  from  the  Hertzian 
vectors  by  setting 

(35)     A  =^^11  -tj/V^l^  -   V^ 

The  associated  scalar  potential  function  is 
(86)     ^  -  -  VJZ  +  H 


d/ 


in  which ^  must  also  satisfy 


(87)  K/'^ij^   -mE   0*  -C> 


Boundary  Conditions  -  Discontinuities  in  the  Field  Vectors 

The  above  equations  have  all  been  based  upon 
the  assumption  of  continuous  and  isotropic  media. 
It  is  now  necessary  to  consider  the  relations  V7hich 
arise  due  to  the  boundaries  between  tv/o  media,  where 
there  occur  abrupt  differences  in  the  inductive 
capacities,^  and  (^,  and  in  the  conductivity,  s   . 

To  aid  in  the  analysis,  it  is  useful  to  postulate 
a  very  thin  transition  layer,  S,  betv;een  the  media. 
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wherein  the  parameters  >^ ,  (?  ,  and<5,  vary  rapidly  but 
continuously  from  medium  1  to  nedlum  2.  Within  this 
layer,  S,  the  field  vectors  and  their  first  derivatives 
are  continuous,  bounded  functions  of  time  and  space. 
Consider  a  small  right  cylinder,  Figure  1,  in  the 
layer  such  that  its  end  surfaces  lie  in  the  surfaces 
of  the  layer  and  separated  by  a  length  A  k,  . 

From  the  volume  integral  of  eouation  (5)  and 

IQ 
the  divergence  theorem  of  vector  analysis  ,  it 

follov/s  that 


(88)  fvB  d/ir  -   fs-n  dcj    =0 


where  the  second  integral  is  a  surface  integral 

over  the  walls  and  the  ends  of  the  cylinder,  /7  being 

the  unit  vector  normal  to    S    ,      If  the  elemental 

area  of  the  end  of  the  cylinder  is  sufficiently  small, 

B  may  be  considered  to  have  a  constant  va.lue  over 

the  end.  Neglecting  differentials  of  higher  order, 

(83)  may  be  approximated  by 

(89)  (  B-n,  +  d-Hx)  ^^  +  contribution  of  walls  =  0 


Fi^    I 


surface  S 


0)   £,   M,    Si 


Not/val  Bounc/arcj  ConGl/t/cn 
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where  (^^      Is  unit  vector  normal  to  the  surface  in 
the  sense  Invard  to  medium  1,  and  tl;,    the  unit 
vector  normal  in  the  sense  inward  to  medium  2. 
Since  the  contribution  of  the  ^.zalls  is  proportional 
to  Ai  »  it  vjill  become  vani shingly  small  v.-hen  the 
limit  ic.  taVien  l\S. -^  0    ,  and  the  transition  layer 
shrinks  to  the  surface  S.   Designating  the  value 
of  magnetic  induction  at  the  point  on  S  in  medium  1 
as  B. ,  and  in  medium  2  as  Bp,  and  pasBins  to  the 
limit  AH^  0      and  Aa -^  0 


(90)       (B^-  BJ  n   =  0 


Thus  the  transition  of  the  normal  corapbnent  of  B 
across  any  surface  of  discontinuity  in  the  medium 
is  continuous. 

Applying  the  same  process  of  reasoning  to  the 
electric  displacement  vector  D,  eouation  (lO)  is 
integrated  over  a  volume  and  its  equivalent  from 
the  divergence  theorem  substituted.   In  this  case, 
however,  since  V-  D  is  equal  to  the  charge  density^  f 
when  it  is  integrated  over 'a  volume  it  becomes 
equal  to  the  total  charge  within  the  volume,  q,  thus 


(91)  j^V-Odnr   -  lotidc^    =^ 
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This  charge, q,  is  distributed  throughout  the 
transition  layer  v.'ith  a  density  p  ,    and  as  the 
cylinder  is  redu.ced  tov/ard.  the  limit,  the  q  must 
remain  the  same  because  of  the  conservation  of 
charge.   Thus 

(92)     (9    '     f"  aP  AO 

As  A  Jl  -*  0  ,    f   must  then"»<»o.   It  is  then  convenient 
to  replace  the  product,  f  a^  ,  by  a  surface  density 
or  charge  per  unit  a.rea.fUJ   .   The  transition  of  the 
normal  component  of  the  electric  displacement  vector, 
D,  across  any  surface,  S,  is  then  given  by 


(93)  (0,-D,)  n  =  oj 


The  presence  of  a  layer  of  char,<5e  on  S  results  in 
an  abr'jpt  change  in  the  normal  component  of  D,  the 
araoijnt  of  the  discontinuity  being  equal  to  the 
surface  density  in  coulombs  per  square  meter. 
It  is  also  necessary  to  investigate  the 
behavior  of  the  tangential  components  of  E  and  H  at 
the  boundary.   To  aid  this  analysis,  consider  a 
rectangle  of  length  ^5  and  thickness  aX   »  instead 
of  a  cylinder,  lying  at  the  surfaces  of  the  tv;o 
adjoining  media.  Let  5^      be  the  area  of  this 
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rectangle  and  Cq    "be  its  closed  contour. 

Taking  equation  (l)  and  integratin,'^  over  the 
surface 

(94)  flV^On.da    +/  ^  nda    =0 


where  n^    is  the  unit  normal  vector  in  the  conventional 

positive  sense,  as  shown  in  Figure  2.   By  Stokes 

19 
theorem  ,  this  first  surface  integral  is  equal  to 

a  line  integral  of  the  vector  E  about  a  closed 

contour  enclosing  the  surface.  Thus  (94)  may  be 

expressed  as 

(95)  f  E-d5     -^  J  ^  n^dci  ^0 

Designate  by  y     and  7^  "^1^©  unit  vectors  in  the 
direction  of  circulation  as  shovm  in  Figure  2.   Then, 
neglectinjg  higher  order  differentials,  equation  (95) 
may  be  expressed  by 

(96)  {f  T,  i-ETiiJ^S   t  contribution  from  ends^-^He^SAH 


rig.  z 


(I)  8>  ^,  ^1 


r, 


h 


V~rrA- 


Tan^eniial  Bovfidari)  Cond/fm 
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But  the  contributions  from  the  ends  are  proportional 
to  Ai   which  becomes  vanishingly  small  v;hen  the  limit 
Is  taken.   Let  ^?  be  the  unit  vector  normal  to  S 
as  shovm,  and  T  be  the  unit  tangent  vector  such  that 

(97)  T  =  n,^r\ 

Also,  it  is  a  vector  identity  that 

(98)  /7„  a;  M-  E  =  n^-n  X  E 

Thus,  v,'hen  the  limit  is  taken  A^-^  0       ,  AS -^  0 

(99)  f7.- [«'(£, -Ej  ^iiAf,Al)\  '-0 

Since  n^^  0  and  the  orientation  of  the  rectangle 
is  arbitrary,  the  bracketed  expression  above  must 
be  zero,  and  so 

(100)  n  ^(t^-Z)  =  -bm  ^  a^ 

And  since  we  have  imioosed.  the  conditions  that  the 
field  vectors  and  their  derivatives  be  bounded,  the 
right  Gide  of  the  above  equation  vanishes  when  a^-^0 
Therefore 

(101)  n^fE,-£,)  ^  0 
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and  it  is  seen  that  the  tangential  components  of 
the  vector  S  through  a  surface  of  discontinuity- 
are  continuous. 

The  behavior  of  the  magnetic  intensity  vector, 
H,  at  the  boundary  can  be  deduced  by  the  same 
process.   Inte^^rating  equation  (2)  over  the  surface 
of  the  rectangle  described  above  and  applyin3  Stokes 
theorem  to  the  first  terra,  there  is  obtained 

(102)  In-ds-  jll-n4a    =/ J-n.d^ 

and  It  similarly  follovrs  that 

(103)  n>^(^,'H^    =    bm  (I?  ^rj)A^ 

The  limit  of  -^j  aJI  as  a/->  (9  v;lll  vanish 

because  D  and  its  derivatives  are  bounded.   But  it  is 
Dossible  that  iiit]  -J a  ^      may  not  vanish,  because 

the  current  I  =-  -J-noAS^i  is  squeezed  into  an 
infinitesimal  layer  on  the  surface  as  the  sides 
are  brought  together,  and  -J-^  ^    ,      It  is  convenient 
to  represent  this  surface  current  density  by  K, 
v;hlch  is  defined  as  the  limit  of  t7"a^   as  A^ 
ar>r)roaches  zero  and  J  becom.es  infinite.   Then 
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(104)    m(H,-H,)  =K 


When  the  conductivities  of  the  adjoining  media 
are  finite,  no  surface  current  can  exist,  because 
S  is  bounded  and  the  product  of  sEa^      vanishes 
wlth^j^.   In  the  analysis  of  dielectric  antennas, 
this  condition  exists,  and  so  (104)  can  be  reduced 
to 


(105)   n  /('H,  -k)  =^ 


In  the  foregolns,  expressions  were  developed 
for  the  transition  of  the  normal  components  of  B 
and  D,  equations  (90)  and  (93),  and  for  the 
tangential  components  of  E  and  (H),  equations  (lOl) 
and  (105).  By  using  (15)  and  (l6),  the  normal 
components  of  H  and  E  are  immediately  obtained  for 
the  transition  of  an  electromagnetic  field  from 
one  medium  to  another. 


(106)   hYh,  -  ^J^)  ^0 


(107)   n.(E, -|^e)= 


CJ 
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and  for  the  tangential  components  of  D  and  B 

(108)  r^  x{0,  -  ^  O)  ^0 


(109)  nA(6, -^0,)-O   ^ 

(for  finite^) 


Povntlnp:'3  Theorem  jLn  rep^lons  free  of  Hysteresis 

Poynting's  theorem- gives  the  relation  betv/een 
the  rate  of  chan,3e  of  energy  stored  in  a  field  and 
the  energy  flow.   This  may  be  obtained  as  a  general 
integral  of  the  field  equations  (l),  (2),  (5),  and 
(lO),   Since  E  "^  ^^s  "^"^^   dimensions  of  ^ov/er 
expended  per  unit  volume,  it  is  apparent  that  this 
may  be  obtained  by  a  scalar  multiplication  of  (2)  by 
S.   Thus 

(110)  Z'V^W  -   E-|f  -E'-tJ 

Likev;ise  enuation  (l)  can  be  given  the  units  of 
pov/er  per  unit  volume  by  multiplying  it  by  H,  Thus 


se 


(111)  H-V^E  +H'||  ^0 
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Subtractins  the  above  two   equations 

(112)  H-VA£  -e-VaH+E-It +H-||=-E-^ 

But,  by  vector  Identity,  the  first  tv;o  terms  of 
(112)  are  equal  to  V-f£xH)  and  so  t'cils  equation  may 
be  written  as 

(113)  V(E'^H)  t  r  J  =  -  E-f^  -  h-^ 

Integrating  over  a  volume  V  bounded  by  a  surface 
S,  and  applying  the  divergence  theorem  to  the  first 
term,  this  becomes 

(114)  i(£*h)-ncla   tjfEJd. -'-{(£■  i^H-f^ 

The  terms  on  the  right  above  represent  the  rate  of 
decrease  of  the  electric  and  magnetic  energy  stored 
In  the  volume.   The  second  term  on  the  left  expresses 
the  sum  of  the  losses  in  Joule  heat  and  the  pov.-er 
exToended  by  the  flov7  of  charge  against  the  impressed 
forces.  If  there  are  no  transformations  of 
electromagnetic  energy  into  elastic  energy  of  a 
stressed  medium,  the  balance  is  maintained  by  a  flow 
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of  electromagnetic  energy  across  the  surface  bounding 
the  volume.   This  is  Poynting's  interpretation  of 
the  significance  of  the  surf^.ce  integral  in  (ll4). 
The  reduction  of  electromagnetic  energy  stored  in 
the  volume  is  partly  accounted  for  by  the  Joule 
heat  loss,  partly  by  the  energy  introduced  by  the 
impressed  forces,  and  the  remainder  flows  outward 
across  the  boun:"lng  surface  S,  which  is  represented  by 

(115)  Isnd^     =l(E^H)-ndi^ 

Now  the  Poynting  Vector  Is  defined  by 

(116)  S-E^H  v/atts/meter 


2 


which  may  be  Interpreted  as  the  intensity  of  energy 
flow  at  a  point  in  the  field,  or  in  other  vrords, 
the  energy  per  unit  time  crossin<5  a  unit  area  whose 
normal  is  oriented  in  the  direction  of  the  vector 
E^H  .   It  is  Poynting *s  theorem  that  the  vector 
product  (ll6)  at  any  point  is  a  measure  of  the  rate 
of  energy  flow  per  unit  area  at  that  point. 

Another  useful  form  of  Poynting* s  theorem  is 
the  differential  form,  in  which  the  Joule  heat 
dissipation  is  separated  from  the  ■cower  ex-oended 
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by   the  flovr  of   charge  against   the  impressed  forces. 

V/hlle  E  re-oresents   the   total   electric  field 

i 

intensity,  let  S  be  the  intensity  of  the  impressed 

electromotive  forces.   Then 

(117)      E'i-E.' 

and 

(118)         £.J  ..^    ,^'.J 

Substituting  (llB)  into  (113)  and  replacing  D  by 
its  equal  ^fTand  B  hy^Hf    "t-he  differential  form 
becomes 

(119) 


The  Gom-ilex  Poyntin^  Vector 

Let  h  ^M-tuj    represent  the  density  of  the 
electromagnetic  energy  at  any  instant  and  p  =gJ^-E'--J 
be  the  rjower  expended  per  unit  volume  in  Joule  heat 
and  by  the  flo-.-;  of  charge  against  impressed  forces. 
From  (119)  v:ith  these  substitutions,  the  Poynting 
theorem  from  (113)  can  be  v/ritten 
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(120)    V5  ^  ^  ^^  ^  0 


In  a  stationary  field,  h  is  independent  of  time  and 
the  middle  term  v;ould  be  absent.  Moreover,  v/hen 
the  field  is  a  simple  harmonic  variation  of  time, 
the  time  average  of  the  energy  density  over  a  period 
is  constant  and   -^  =0       ,    the  bar  over  the  symbol 
indicating  an  avera^^e  value  for  a  period.   In  the 
case  of  periodic  fields,  it  follov/s  that 


(121)  or 


5 

Now  these  field  vectors  are  harmonic  functions  of 
time  and  can  be  represented  in  comolex  form  v/ith 
real  and  imaginary  parts.  By  complex  variable 
theory,  it  can  be  shovm  that  the  time  averaf^e  of 
the  product  of  the  real  parts  of  two  complex 
quantities  is  equal  to  one  half  the  real  part  of 
the  -product  of  the  first  quantity  times  the 
conjugate  of  the  second  quantity.   Symbolically 


(122)  (^Mj-^MJ  =i(^(A>A^) 
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where  A  is  the  comolex  conjugate  of  Ag. 

According  to  the  above,  the  mean  Intensity  of 
the  energy  flow  in  a  harmonic  electromagnetic  field 
is  the  real  uart  of  the  complex  vector  2  EXH 


5   =  (^(£)/^a^(H) 

This  is  the  complex  Poynting  vector  v.'hlch  may  also 
be  denoted  by 

(124)     S*  =  iExH 

In  a  plane  v/ave  the  magnitude  of  H  is  -§-  »  v;here 
/j  is  the  impedance  of  the  medium.  Therefore  in 
engineering  practice   the  ex-nression 


(125)  Power    ^    -|  ^ 


is  frequently  used,  corresponding  to  the  vector 
expression  (124). 
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Representatl on  of  a  cylindrical  field  by  Hertz  Vectors 

In  a  homogeneous  medium  v/hich  Is  free  from 
any  fixed  polarization,  the  Max^vjell  Equations  (l) 
and  (2)  are  satisfied  by  the  follouing  relations  in 
terras  of  Hertz  vectors 


(126)    e:  =  v/v^n  -mv^  If' 


The  first  of  these  can  be  readily  verified  by  taking 
equation  (Bo),  applying  the  vector  identity  (38) 
to  the  first  term,  and  noting  that  v^J-^^-^^=^ 

"by  (67).   The  second  one  is  similarly  verified. 
Now  consider  a  set  of  coordinate  surfaces, 
Figure  3,  formed  by  a  family  of  cylinders  whose 
elements  are  parallel  to  the  ^  axis.   These 
cylindrical  surfaces  need  not  necessarily  be 
circular  or  closed.   Unit  vectors  i,  ,  i^  ,  and  1_ 
are  assigned  as  shovm,  i-  being  normal  to  the 
cylinder,  1^  tangent  to  it  and  parallel  to  the 4 axis, 
and  lo  is  tangent  to  the  surface  and  perpendicular 
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PIGURS  3. 


Relation  of  Unit  Vectors  to  a  Cylindrical  Surface, 
Tho  'generators  of  tho  surface  are  parallel  to  the 

£  axis. 
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to  i   and  1  ,  forrain.'3  a  rlc;ht  handed  orthogonal 
set.  Position  v/ith  respect  to  axes  determined  by 

the  three  vectors  is  measured  by  the  coordinates 

1   o 
A   ,  M'    ,     ^      .   An  Infinitesimal  line  element  is 

(128)  dd-U,dM'  ^lJJm"-  fl.dB 

Consider  an  electromafrnetlc  field  in  an  isotropic, 
homogeneous  medliim,  unbounded  in  extent,  existing 
such  that  the  Hertz  vector,  II,  is  directed  along 

the  H  axis.   11^  =0   H^  =  0,  and  11^  -  (/'  . 

■Ji- 
ll =0.   The  electric  field  intensity  vector  of 

this  field,  as  ^iven  by  (126)  is 

(129)      E^'^  =  vxvxjr 

and  the  magnetic  vector  of  the  field,  by  (127)  is 


(130)  H^'^^CSh-  ^^)vxjr 


The  components  of  the  vectors  E  can  be  calculated 
from  (129)  by  means  of  the  vector  formula  for  the 
curl  in  an  orthogonal  system  of  curvilinear 
coordinates  (reference  19,  page  90  ),  applying  this 
tvrice. 


46. 


(131) 

^hen     A'-%        -^.-^         ^"% 


Since  II,  =  JJ^^  =  (9   ,  this  comes  out  to  be,  after 
takin*'^  the  curl  twice 


(132)  £r.^  =  4i!^ 


Likewise,  to  compute  the  components  of  H  from  (I30), 
simply  take  the  curl  once  and  multiply  by  £[^  '^  ^^  ^ 
thus 


1 


^ 
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(133)   ^f-(4-Uff 


From  the  z  scalar  comnonent  of  a  function 
U^  - y^ ,  an  electromagnetic  field  has  been  developed 
which  is  characterized  by  the  absence  of  any  axial 
or  longitudinal  component  of  the  raa.^netlc  vector. 
It  vjas  earlier  Dointed  out  that  II  is  the  electric 
polarization  potential,  and  therefore  this  field 

may  appropriately  be  called  a  field  of  the  electric 

A 

type,  but  follovring  the  terminolorjy  of  Schelkunoff  , 

it  is  more  often  referred  to  as  the  "transverse 
magnetic  field," 

It  is  necessary  to  digress  to  develop  here 
the  wave  equation,  ;\^hich  the  function  E^  =  ^ 
must  satisfy. 

Taking  the  curl  of  equation  (l)  and  substituting 
mH     f  or  B 

(134)    Vx(7X£-;  =  -^^(^X^) 
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Since  J  =  6"iE"  and  D  =  £  E"  ,  equation  (2)  can  be 
written 

(135)  VxH=<5E-^^-^ 

Substitute  (135)  into  (134)  and  making  use  of  the 
vector  identity  (AO) 

(136)  v(v-£)  -  v-vr  =  -^{^0^"^^|| 

Rev7ritin^  equation  (lO),    substituting  ft   for  D 

(137)  V'E    -  -|- 

Since  in  free  space  or  in  a.ny  conducting  medium  f 
is  independent  of  the  field  distribution,  the  first 
term  of  (136)  may  be  taken  to  be  equal  to  zero. 
Then  (136)  becomes 

(138)  ^%-^<^0  --^^^  ^0 

Since  ^=^£   ,  this  could  also  be  written 
(139)    r7'j-^s^.--^<S^   =0 
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By  elimlnatini3  E  from  Max\;ell*s  equations,    there  is 
also  obtained 


(140)       vV._^^|^   _^g||^, 


Now,  returning  to  the  equations  of  the  cylindrical 
field,j^  =  ^,  the  components  of  E  and  H  ^iven  in 
(132)  and  (133)  must  satisfy  the  scalar  v;ave  enuation 


(141)  vV  -^cf  0  --^^^  -^ 


The  Laplacian  of  the  scalar^,  expressed  in  K;eneral 
orthogonal  coordinates  (reference  19,  page  ISO     ),is 
(142) 

Substituting  this  in  (l4l),  \ihere  A-.   =  1,  nf=i 
(143) 

The  elementary  harmonic  solutions  of  (l43) 
are  of  the  form 


t 
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/  ,,>)     *^j<i-^'*'^ 


(144)        I//  =  {(u'^u'')  e 

vrhere     -f(u'jU^)  is  a   solution  of 

(145) 


Similarly,  a  partial  field  can  be  derived  from  a 
second  Hertz  vector  II  by  the  operations 


(146)   5^'^--^fc  V^l"" 

If  the  vector  II  is  directed  along  the  2  axis, 
the  vector  components  come  out  to  be 
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til  ^  9i  d<^ 


(^  _  ±    a'jt 


(1^8)   Hf'-±M^ 


*x1 


The  scalar  function  II^  is  a  solution  of  (l4l), 
and  is  a  field  of  the  magnetic  type  or  transverse 
electric,  characterized  by  a  zero  longitudinal 
component  of  E, 

If  the  partial  fields  derived  from  II^  and  11^ 
are  superposed,  the  electroraac'^netic  fields  are  of 
sucljL  generality  that  they  can  satisfy  any  prescribed 
set  of  boundary  conditions  on  any  cylindrical  surface 
vhose  generating  elements  are  parallel  to  the  Z 
axis.   In  general,  the  choice  of  the  coordinates  is 
limited  practically  to  systems  in  which  the  variables 
of  equation  (l45)  are  separable.   It  is  for  this 
reason  that  this  development  of  equations  of  electro- 
magnetic waves  in  a  rod  is  limited  to  the  case  of 
the  circular  cylinder. 
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V/ave  Functions  of  the  Circular  Cylinder 

The  simplest  of  the  separable  cases  is  that  in 
v;hlch  the  family  (^'=constant  is  represented  by  a  set 
of  coaxial  cylinders,  v;here  M'=r  »  U^--d      $    ^,  - 1 
and  -^i-f.      See  Figure  4.   In  this  case  equation 
(l45)  reduces  to: 
(149) 

v;hich  can  be  separated  by  v/riting  f(l^,9)^^  ^   product 

(150)      f  =-f:fr)  fje) 

vrhere  f/r)is  an  arbitrary  solution  of  the  eouation 
(151) 

and  fa (^j  is  an  arbitrary  solution  of 
(152)      ^   +^V,  =  0 

Solutions  of  (151)  are  Bessel's  functions,  and 
a  solution  to  (152)  is  f>(6)  -  ^^^^  ,  The  nararaeter,  n, 
is  a  separation  constant.   Its  choice  is  governed  by 
the  physical  rsquireinent  that  the  field  must  be 
single  valued  at  any  fixed  point  in  space.   Since  at 
the  present,  the  discussion  is  limited  to  continuous 
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Ax 


Coordinates  of  the  Circular  Cylinder 
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media,  the  field  is  periodic  in  9  and  the  values  of 

n  are  limited  to  integers,  0,  1,  2,  3 ,  either 

positive  or  negative,  Hov/ever,  in  a  field  represented 
by     -oarticular  solutions  of  (151)  and  (152),  which 
is  bounded  by  planes  O^O^  ,    and  (^  =  (9j^  ,  non-integral 
values  V7ould  in  general  have  to  be  assigned  to  n. 

The  constant  h  is,  in  ,c;eneral,  a  com-clex  quantity 
pertainln,':^  to  the  propaj7;ation  of  the  uave  and  is 
called  the  propagation  factor.  An  explicit  expression 
for  h  in  terms  of  the  frequency  and  the  constants 
of  the  medium  can  be  obtained  only  after  "he  behavior 
of  (//    over  a  cylinder   T^  constant  on  a  plane  ?= 
constant  has  been  prescribed. 

The  equation  (151)  Is  Bessel's  equation,  the 
solutions  of  •'.rhlch  are  the  Bessel  functions  -rLn  (if^'^ r) » 
In   accordance  with  terminology  emx-)loyed  by  Jahnke 
and  Emde  and  several  other  authors  the  term  "circular 
cylinder  fionction"  v;ill  be  used  to  denote  any 
particular  solution  of  (151)  ^.nd  designated  by  the 
letters  f|  -  Z^  (IS*^  r)   .  The  order  of  the  functions 
whose  argument  is  ^^^^  f     is  n.  Particular  solutions 
of  the  scalar  v:ave  equation  (l4l)  v/hich  are  periodic 
in  t  and  6  can  be  constriacted  from  elementary  waves 
of  the  form 
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(153) 

Reference  is  made  to  standard  texts  in 
mathematics  for  the  properties  of  the  various  Bessel 
functions,  but  it  is  essential  to  this  development 
to  note  several  pertinent  characteristics.   The 
Bessel  function,  or  circular  cylinder  function  of 
the  first  kind,  is  finite  and  converges  for  all 
finite  values  of  the  arj^ument.   For  non-integral 
values  of  n  it  has  tv'o  linearly  inder^endent 
solutions  for^(f)        and  -xL^(f)      ;  but  if  n  is  an 
integer  the  tv/o  solutions  are  not  Independent, 
A  second  inde^^endent  solution  can  be  obtained  by 
the  Bessel  function  of  the  second  kind,  /\,i^  (^)        , 
defined  by 

( 154)  U„ (!)  =  jbr  C^-.  (f)c^'"Tr  -^(f)] 

When  n  is  an  integer,  N^(f)  Is  a  second  solution, 

but  this  is  useless  at  f~0,    for  the  function  becomes 

infinite  at  this  point.   Solutions  can  be  obtained 
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in,  terms  of   still   other  functions,    ttie  Hanliel 
fiinctlons,     H^  (f)        and      ff^  (f)       ,    v/hich  are  merely 
linear  combinations   of   the  first  and   second  Bessel 
functions  as  follows 


(155)  Hmif)  =^^(f)+ANjf) 


^    (156)      Hi^(f)  =^J?)-,HJ?) 

But  the  Hankel  functions  likewise  become  infinite 
for  f-0  and  so  are  excluded  from  use  where  the 
argument  is  very  small.   Expansion  of  the  J  and  N 
functions  shov/s  a  similarity  to  the  cosine  and  sine 
functions  and,  as  such,  represent  standin,-;  waves 
when  used  with  the  time  factor  e''^^      .      The  H 
functions  correspond  to  exponential  functions  v.dth 
imar^inary  exponents,  and,  when  used  with  the  time 
factor  e''^^'^       represent  traveling  cylindrical  \ravea, 


The  Field  of  Circularlr  Cylindrical  V/ave  Functions 

Every  electromagnetic  field  v/ithin  a  homogeneous 
Isotropic  domain  can  be  represented  by  -linear 
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combinations  of  the  follov/lng  functions 


^^57)  f^^^=e'^'^^(W^r)e'"'^-'"* 


^^^^^  U^  =  e^"' H^'  (i^^^*-^  e-^^-'"'^ 


The  first  equation,  (157),  alone,  is  sufficient  to 
represent  an  electromagnetic  field  v/ithln  finite 
domains  including  the  axis  fr^  .  But  at  great 
distances  from  the  source,  (158)  must  be  used,  for 
B.Q\f^^^T)       becomes  very  large,  the  function  reduces 
to  a  wave  travelling  radially  outv/ard.   This  can  be 
more  readily  understood  v/hen  it  is  noted  that  for 
very  large  arguments  (f)    the  Hankel  function  is 
approximately  given  by 


Each  elementary  wave  is  identified  by  the  parameters 
n,  h,  k.  V/hen  n  is  zero  the  field  is  independent  of  6 


I 
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and  Is  synunetrlc  about  the  8.xls,  v/hen  h  is  zero  the 
propap;ation  Is  purely  radial  and  the  field  becomes 
only  tvo  dlmenv^-ional.   The  tw.o  functions,  (157)  and 
(158)  may  represent  inhoraof;eneou3  plane  waves  as 
v/ell  as  homogeneous.   The  planes  of  constant  phase 
are  propa.^ated  along  the  ^  axis  with  a  velocity 
\/-^  V7here  cA  Is  the  real  part  of  h,  but  the 

amplitudes  over  these  Planes  are  functions  of  f 
and  0. 

The  formulae  for  the  components  of  the  field 
vectors  can  be  obtained  from  equations  (147)  ^-nd  (l48), 
remembering  that  II  =11^-0  and  Ilg  =  ^      , 
Also,  the  above  equations  were  for  a  generalized 
cylindrical  field,  while  we  are  now  restricted  to 
circularly  cylindrical  functions  and  therefore,  A^-Atrlt 
J^;^^f     ,      and  ^'=r     )<A^-d    ^  ^  ~Z  •  Note  also 

that,  due  to  the  exponential  ^ -^■^t-'<-(*^i       ,  differen- 
tiation v/ith  respect  to  2  multiplies  the  function 
"by   1^'A   •  ^or  ^^'^   transverse  magnetic  field,  the 
results  are! 
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(161) 


H^ 

.^co   r  de 

H«» 

Hf  = 

0 

similarly  for  the  transverse  electric  field 


(162)  E':''-.,^^r 


It-  0 


n^    -  -aJL  ;)  \r 


dr 


(165)   H«^-*^!f 


ht^  =  (r-A')f 
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V'hen  initial  conditions  are  prescribed  over  given 
plane  or  cylindrical  surfaces,  a  solution  can  be 
constructed  by  suDerposing  these  elementary  wave 
'functions.   For  fixed  values  of  W  (okJ^)  and  h, 

the  field  in  cylindrical  coordinates  may  be 
represented  by  the  equations  (164)  and  (155)  below, 
in  which  the  d^-^   and  ^  coefficients  must  be  determined 
from  the  initial  conditions. 


(164)  £,  = -^  X  .<^^. -.^.  f  fe,f 


/>?=-  oo 


^^  =  i?  5^^^--^^$^ 


dm 


(165)   He-£  2  <?.^^-4t-'i.^ 
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Pt^opan^atlon  aloiy>^  a  circular  cylinder 

Consider  a  circular  cylinder  of  Infinite  length, 
of  radius  a,  of  material  v/ith  propagation  constant 
kl  imbedded  in  an  infinite  homogeneous  medium  of 
proDai^ation  constant  k^  ,   The  field  muct  be  finite 
at  the  center  and  consequently  the  v;ave  functions 
vjithin  the  cylinder  will  be  constructed  of  Besael 
functions  of  the  first  kind.   Outside  the  cylinder 
the  field  will  extend  to  Infinity  and  consequently 
the  functions  must  be  constructed  of  Hankel  functions. 
In  the  equations  below  X/*ls  substituted  for  (^f-A^) 
and  xl    is  substituted  for  (^/-X)   .   The  function 
^*^e -t^^2  ^Auj;t  is  abbreviated  to  Fn.   Summations 

are  in  each  case  from  n=i-°^ton=»t«3      •  A 
prime  above  a  cylinder  function  indicates  differentiation 
with  respect  to  the  argujnent  [X  0   .   The  super- 
scripts 1  or  e  over  the  coefficients  ^^       and  b^ 
indicate  points  interior  and  exterior  to  the  cylinder 
respectively. 


■=  igj^f-f-J^i-^^M^ij 


r, 


'K 


(166)      £j^-2   [ff,^(^/)«.H  '-f^Xp(:')l>i]f^ 


r;=  2  &'Maar^ 
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H>  2  [^,^.rxo^;+^4:fx,r)^^]r^ 


(167)      K=   ^l^x^^M'^^^-'^M'O^^ 


r, 


nn 


h;=  2  p-':v-;A;]r.. 


e;=  2  [f  <^V^o^.^-i^^Hi'^M^f 


u 


(168)  e:=-2  r#,Hiv^^)^.^^^-f^/iyMdF^ 


xi 


E'  =  2  LH^''M^lJ/5^ 


HN2  [^,H^'^M^.^^f/^/(-^^O^f 


fA 


M 


(169)      H!=   5 


[J^rty^.}.-f,H:'((jrj«;; 


i/ti. 


H;=  ICH^'MV^JF^ 
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The  coefficientB  of  the  expansions  and  the 
propaf^atton  factor  h  muct  next  be  determined.   At 
the  bo'ondary  v;here  r  is  equal  to  the  constant  radius 
of  the  cylinder,  a,  the  tangential  components  of  the 
field  are  continuous,  v;hich  condition  supplies  a 
relation  betv;een  the  coefficients.   From  the 
continuity  of  the  tangential  coniDonents  of  E,  the 
follovring  relations  are  found,  where  the  substitutions 
have  been -made  X,a=P  ,  Xi^-AT 

(170)  ' 


(171)   J„(f)a:;    =  \ti^)a.' 

From  the  tansential  components  of  H,  similar 
relations  are  found  as  follov/st 

(172) 
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(I) 


e 


(173)    -J.Mb^     =  H;fu}4 

The  above  relations  constitute  a  honof;;eneous 
system  of  linear  equations  satisfied  by  the  four 
coefficients  c?^  ^  bJ:  ,  ^J  ^  b^  .   There  vill  be 

a  nontrivial  solution  only  if  the  rank  of  the  matrix 
Is  less  than  four,  vjhich  i-'ill  be  true  if  its 
determinant  is  zero.   The  propagation  factor  is 
determined  by  the  condition  that  the  detorminant 
of  (170)  throu-h  (173)  shall  vanish. 

Expansion  of  this  determinant  leads  to  the 
transcendental  equation 

(174) 


The  roots  of  (174)  are  the  allov/ed  values  of  the 
propaf^ation  factor,  h,  and  so  determine  the 
characteristic  or  natural  modes  of  oropa.^ation. 
Thus,  there  has  been  developed  an  eouation 
v.'hose  solution  enables  calculation  of  natural  modes 
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and  the  propao;atlon  constants  for  long  rods.   Equation 
(l49)  expresses  the  true  result  of  Maxi^ell  in  aji 
infinitely  loa^  rod  imbedded. in  a  pei-^ectly 
homogeneous  medium.   Equations  (157)  and  (158)  are 
valid  solutions  to  this  equation.   By  use  of  the 
continuity  of  the  tangential  components  of  E  and  H 
at  the  boundary  of  the  rod,  it  v/as  possible  to  use 
(157)  and  (158)  to  find  relations  which  enable  one 
to  solve  simultaneously  for  the  propagation  constants. 

Solutions  to  (174)  have  been  obtained  and  it 

9  10 
has  been  shown  *   that  they  verify  exoerimental 

results  v/ithln  the  accuracy  of  experimental  procedure. 

In  the  next  chapter,  a  nomogram  will  be  presented 

which  provides  solutions  for  the  ca^e  of  dielectric 

rods  in  air  or  vacuum. 

It  should  be  emphasized  that  the  preceding 

development  assumes  that  the  rod  is  long  enough 

to  neglect  boundary  conditions  at  the  ends,  and  that 

the  radius,  a,  is  not  microscopic  as  far  as  the 

external  v^ave  is  conceirned.   The  question  as  to  v^hat 

size  rod  is  required  to  Justify  the  assumption  that 

the  radius  is  not  microscopic  is  most  pertinent. 

If  the  radius  is  small  compared  to  a  v/avelength, 

it  is  probably  in  error  to  assume  that  Hankel  functions 
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are  required  for  all  v/aves  exterior  to  the  rod 
while  Bessel  functions  of  the  first  kind  suffice 
for  internal  waves. 

Havins  obtained  the  propaj^ation  constants, 
available  published  literature  indicates  that  all 
authors  proceed  to  determine  radiation  characteristics 
on  the  assuinr>tion  that  the  propagation  is  truly 
constant  throa^hout  the  rod.   The  experimental 
evidence  that  this  is  not  true  in  short  rods  of 
diameters  of  a  wavelength  or  less  has  apparently 
been  neglected,  and  available  published  theories 
on  pov/ers  show  serious  defects  when  confronted  with 
experimental  evidence.   For  examiole,  experimental 
curves"*    shov;  that  Dov/er  on  axis  is  linear  with 
exposed  length  of  dielectric  for  short  antennas. 
The  published  theories  do  not  permit  such  a 
relation, 

A  satisfactory  solution  to  this  problem  from 
a  theoretical  standpoint  cannot  be  completed  in 
the  present  v:ork.   The  experimental  data  presented 
may  provide  a  stepping  stone  toward  its  ultimate 
achievement. 
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CHAPTER  2. 

APPLICATION  OF  HXECTROI-IAGNSTIC  TK20RY  TO 
CIRCULAR  DI^TLECTRIG  RODS 


In  Chapter  1  an  9quatlon  was  developed 
expressing  In  circular  Coordinates  the  nature.l 
modes  of  electromartnetlc  v;aves  in  an  infinitely  long 
cylinder  imbedded  in  an  infinite  homogeneous  medium. 
The  next  step  is  to  specialize  this  to  the  case  of 
a  dielectric  cylinder  imbedded  in  air  or  vacuum. 
Furthermore,  particular  solutions  of  (l?^)  are 
desired,  and  to  achieve  this,  the  equations  will  be 
put  in  a  form  more  suitable  for  use  with  available 
tables  of  cylinder  functions. 

The  equation  (17^)  niay  be  expressed  in  Modified 
Bessel  functions  (K   functions)  by  use  of  the 
relation 

2i  , iC) 


(175)      Kjr)=xfe"*'  C(^f) 
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or 


and 


(m)      H'fM-fe-'^CU) 


Letting  -iv  rs.   q,  then 


V/e  let  the  above  expression  be  -g,  by  definition, 
to  simplify  notation.   Also  let 

Substituting  these  relations  in  (l?^) 
(180) 
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Since  we  are  only  concerned  with  dielectric  rod 
and  air  or  free  space,  \-:e   may  consider 

(181)   x<^^  M:,    -^  M. 

The  constants  k  and  h  vrere  introduced  in  our 
elementary  harmonic  solutions  of  the  scalar  wave 
equation  in  a  cylinder,  equations  (143-145).  From 
this,  k  v;as  equal  "^(x^^m. S ->i uj6m.  ,  and  we  may  now 

set  S-0         ,    since  the  conductivity  of  dielectric 
and  air  are  negligible.   Thus  JJ^  -  ^UiT,  ,and 

-^3L~  I^'^mJI       ,  In  the  notation  introduced  v;ith 
equations  (166-169),  xf -^^--A""  and  X^^A^-A''    , 
Since  X,=  -^  and   X^"=?=^,  it  follows  that 


(182) 


(1B3) 


V/here   the   subscript  on  £^  has  now  been  changed  to 
zero  to  indicate  the  s"oecial  case  of  free  sr^ace. 
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Eliminating  h   from  the  above,  the  following 
is  obtained 


(184)       -^—f"^  -  (Aj'M.{i,-io) 


Let  us  now  set  for  the  relative  dielectric  constant 
(185)     f 


and  use  the  relation 


(186)     C  ^  ^ 
to  put  equation  (184)  in  the  form 


(187)    T    -    -17     -^TT- 


V7hich  v;ill  be  one  of  our  principal  relationships 

to  use  in  solvin-^  (174)  for  the  propagation  constants. 

To  obtain  the  other  principal  equation,  v;e 
substitute  into  (l80)  the  value  of  h  in  terms  of 
the  other  constants  as  -^.iven  in  ecuatlon  (182), 
substitute  for  k^  and  k  their  corresponding  values 
as  above,  and  substitute  for  u)   from.  (187).  After 
simplification,  using  (186),  the  result  is 
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For  a  given  n,  u)  ,    and  £,  the  last  two  ecuatlons 
can  be  solved  for  p  a»d  q  ^vhich  define  the  mode 
propagated  alon^  the  rod,  and  enable  calculation 
of  the  other  constants. 

From  equations  (l82)  and  (I83)  some  interestiixs 
deductions  can  be  drav/n  concerning  the  longitudinal 
propagation  factor,  h.  V'e  note  from  (I83)  that  if 
q  is  zero,  that  his  equal  to  cuYm^o    »  and  therefore 
the  v/ave  is  traveling  vfith  the  velocity  and  wavelength 
characteristic  of  free  space.   If  p  is  zero,  (182  ) 
shoves  that  h  is  enual  ooS^uJ]       ,  and  the  v;ave  is 
traveling  v;ith  the  wavelen,gth  characteristic  of  the 
dielectric.   For  a  given  dielectric,  freauency, 
and  diameter,  eauation  (I87)   shov/s  that  the  sum 
(p  +  Q  ^)  is  a  constant.   It  will  later  be  shovm 
that  p  is  subject  to  less  variation  than  q.   If 
q  is  small,  which  Is  the  case  for  small  diameters, 
the  K  (q)  becomes  very  large  compared  to  J  (p). 

Since  the  wave  functions  outside  the  cylinder  v;ere 

(1) 

constructed  from  H     functions,  which  v/ere  later 

n  ' 

converted  to  K^^  functions,  it  is  apparent  that  the 
enei»gy  outside  the  dielectric  is  large  compared  to 
that  inside.   Also  for  small  q,  the  vravelength  and 
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velocity  characteristics  are  near  that  of  free 
sDace.   Conversely,  for  large  diameters,  q  is 
large  relative  to  p,  and  K^(q)  becomes  very  small 
compared  to  J  (p) •   Thus  the  energy  is  contained 
mostly  Inside  the  dielectric,  with  a  relatively 
thin  film  of  energy  on  the  outside  surfaces.   The 
v/avelength  in  this  case  approaches  the  characteristic 
of  the  dielectric.   Qualitatively,  vie   mi^^ht  expect 
In  vievf  of  this,  unless  the  dielectric  is  ideal 
and  lossless,  that  losses  Increase  rapidly  v.'ith 
Increasing;  diam.eter  of  the  rod« 

An  important  factor  needed  for  the  design 

of  dielectric  antennas  is  the  ratio  of  free  space 

of 
wavelength  to  the  wavelength 'ipropan;ation,  or 

appa,rent  index  of  refraction.   V'e  have  noted, 

where  conductivity  is  zero,  that! 

(a)  the  propagation  constant, characteristic 
of  the  dielectric,  k,  ,  is  equal  cj YmSi 

(b)  the  propagation  constant  characteristic 
of  free  space, 

(c)  the  l.pngitudlnal  propagation  factor  for 
the  combination  is 

h  =  -^ 
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It  follows  that 


(189)     ^^  ~  1^  -  ^' 


where  n„   is  used  as  a  symbol  for  the  apparent 
a 

index  of  refraction.   In  order  to  obtain  an  explicit 

expression  for  n   In  terms  of  the  other  constants 

a 

in  these  equations,    divide    -jSquation  (130)   by  k 


and  note  that 


Solving  for  ^T.   yields 


xi~ 


(190)     A^  =.  ^^  ^  (e^JiXf-^^i 


Substituting  for  the  numerator  its  quivalent  given 
by  (188),  and  extracting  the  square  root,  we  have 


n   in  terms  of  €  ,  p,  and  q» 


'/>. 


(191)      ^  =  p^rV,^ 
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uhlch  is  In  a  convenient  form  for  obtaining  n 
vjhen  p  and  q  have  bean  determined. 

Curves  of  calculated  a-oparent  indices  of 
refraction  for  various  diameters  and  dielectric 
constants  for  the  dominant  mode  (n  =  1)  are  given 

in  Fi,';-ure  5.   These  are  in  reasonably  close  agreement 

1  p 
with  experimental  values  obtained  by  McKinney  . 

It  can  be  readily  understood  that  this 

a-oparent  index  of  refraction  is  of  fundamental 

importance  in  the  design  of  a  dielectric  a.ntenna. 

It  has  been  shovm  '   *    that  the  optimum  length 

of  a  cylindrical  dielectric  rod,  to  maximize 

radiation  on  axis,  is  that  length  v.'hich  causes  the 

wave  to  emerge  from  the  end  delayed  by  a  half  period 

behind  the  external  vrave.   This  length  has  been 

s hovrn° » ^ ^ » ^3  to  be 

U92)     Jl,    =  i^ 

which  is  obtained  as  follows.  Consider  that  the 
time  for  the  v;ave  phase  to  pass  through  the  dielectric, 
t,,  must  be  greater  than  the  time  for  the  phase  to 
sv.'eep  through  the  same  distance  in  air,  tg,,  by  one 
half  period,  7^  . 
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Substituting  c=  -^X  and  v^  =  ^^    =  liL  2^ 

and  solving  f or  ^o »  ©ouatlon  (192)  is  obtained. 
The  above  expression  for  the  optimum  length  has 
been  verified  experimentally  by  several  investi- 
gators'^* ■'■•^. 

In  the  discussion  so  far,  the  mode  of  the 
electromar;netic  v/aves,  v/hich  is  determined  by  the 
n  in  the  Drecedin^  equations,  has  not  been 
specified  and  the  equations  hold  for  any  mode. 
It  V7as  sho^-m  in  Chapter  1  that,  if  the  wave  is  to 
be  periodic  in  9,  that  n  must  be  an  integer*.  '  For 
a  practical  dielectric  antenna,  it  is  necessary 
to  specify  this  further.   The  node  normally  used 
for  propagation  in  the  wave  guide  and  for  exciting 
the  dielectric  is  the  so  called  dominant  mode, 

V7hich  is  the  mode  of  propagation  having  the  lovmst 

4 

cut  off  frequency.   Shelkunoff  has  sho^-.Ti  ,  by 

allov/ing  q  to  approach  zero  and  evaluating  the 
remaining  terras  In  the  equation,  that  the  n  =  1 
mode  has  a  cut  off  frequency  which  is  theoretically 
zero.   For  n  zero  or  greater  than  1,  it  is  also 
shoi^m  that  the  cut  off  is  not  zero.   Therefore, 
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the  dominant  mode  is  n  *  1  and  thiB  is  the  only 
mode  considered  for  the  remainder  of  this  work. 

The  method  of  solution  of  equations  (l-?) 
and  (188)  must  be  either  numerical  or  graphical, 
and  "both  of  these  methods  are  cumbersome  and 
tedious.   In  order  to  facilitate  obtaining  ra.-:>id 
solutions,  at  a  small  sacrifice  of  precision, 
nomogra:oh  c'.iarts  A  and  B  have  been  ":!repared.   See 
Figures  6  and  7. 
Chart  A,  Figure  6,  is  a  nomogram  of  the  equation 

P    2 

p  +  q  =  constant 

and  Chart  B,  Figure  7,  is  a  nomogra.m  of 

To  obtain  a  solution,  \ie   seek  values  of  p  and  q 
to  satinfy  both  equations.   It  is  assumed  that  the 
relative  dielectric  constant,  6 ,  the  radius,  a,  and 
the  free  space  wavelength,  /[o   ,  are  knovm.   To  use 
the  accompanying  charts,  proceed  as  follows: 

Calculate  the  constant,  p'-  f  q^  ,  by  slide  rule, 
from  the  equation  (187).   Thus 
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and  locate  this  point  on  the  center  column  of 
Chart  A.   A  straight  ed^e  nivotod  at  this  point 
will  give  the  values  of  p  and  q  allowed  for  the 
given  dielectric  and  diameter.  Assume  a  value  for 
q  and  place  the  straii^ht  edc^e  at  this  point  on  the 
left  column,  keepin.-^  it  also  at  the  pivot  point 
of  the  center  column.   The  corresponding  value  of 
p  is  read  vrhere  the  straight  edge  intersects  the 
right  column. 

Now  look  at  Chart  B,  and  mark  off  or  visualize 
an  C   line,  corresponding  to  the  given  value,  running 
straight  and  parallel  to  the  dash  dot  lines.   Find 
the  assumed  value  for  q  along  this  line,  as  indicated 
by  the  dashed  q  lines.  At  this  -noint,  thus  located, 
read  the  corresponding  value  of  p  by  eye  inter- 
polation betv;een  the  solid  p  lines.   Compare  the 
values  of  p  as  given  by  Chart  B  with  that  obtained 
from  Chart  A.   If  the  value  of  p  a.~  given  by  Chart  B 
exceeds  that  given  by  Chart  A,  then  the  assumed 
value  of  q  vras  too  large.  Assume  a  lower  value  of 
q  and  again  find  the  corresponding  p  values  from 
Charts  A  and  B  as  before.   w'hen  the  two  charts 
give  identical  p*s  for  the  same  q,  the  solution 
is  correct. 
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The  above  Is  a  trial  and  error  method,  but 
it  is  easily  r^erformed  and  converges  very  quiclcly 
to  a  solution  "".jhich  i3  sufficiently  accurate  for 
many  Dur-ooses.   It  may  liappen  that  a  solution  nay 
be  required  that  does  not  fall  on  the  chart.   It 
was  im^->os5ible  to  make  a  chart  i;hich  was  readable 
and  of  reasonable  size  without  considerable 
restriction  on  the  range,  Hov/ever,  the  range 
selected  for  the  accorar^anylng  charts  is  believed 
to  be  useable  for  most  of  the  practical  problems 
anticipated  at  this  time.   If  the  required  solution 
does  not  fall  on  the  ch-arts,  numerical  methods  should 
be  used.   In  order  to  facilitate  these,  the  values 
of  f  and  g  as  functions  of  p  and  q,  respectively, 
may  be  obtained  from  the  plots  of  these  f^unctions, 
given  in  Figures  8  and  9.   These  plots  may  also  be 
useful  to  refine  the  gra-nhic  solutions  obtained 
from  the  Charts  A  and  B,  if  additional  accuracy  is 
desired. 

Having  solved  the  enuations  for  the  v.'aves  in 
a  circular  dielectric  rod  in  air,  we  can  novj 
readily  calculate  the  optimum  length  for  a  given 
dielectric  constaait,  diameter  ,  and  freouency 
(equation  192),   The  next  problem  is  to  determine 
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what  gain  can  be  expected  from  the  antenna.   There 

is  no  acce-Dted  method  for  calculating  the  gain, 

11  15 
althouj-^h  several  theories  have  been  r)roposed  *    . 

It  is  for  this  reason  that  v/e  turn  to  experiment 
to  obtain  the  gain  to  be  expected  for  the  given 
conditions.   The  next  chapter  v/ill  describe  the 
equipment  and  procedure  used  to  obtain  the  needed 
experimental  data-.  After  describing  these 
measurements,  the  data  v/ill  be  discussed  and  an 
attempt  made  to  correlate  the  results  vrith  Marsten's 
theory  . 

It  will  be  noted  that  the  preceding  theory 
has  dealt  exclusively  v/ith  the  circular  cross  section 
rod  antennas.   The  case  of  the  rectangular  cross 
section  is  much  more  difficult  due  to  the  fact  that 
there  are  discontinuities  on  the  edges,  which  are 
cumbersome  if  not  impossible  to  account  for.   Further- 
more, there  is  no  convenient  coordinate  system  in 
v/hich  to  set  doim  the  mathematical  expressions. 
No  published  theory  of  electromagnetic  vraves  in 
'dielectric  rods  of  rectangular  cross  section  is 
available,  and  consequently  for  the  de.sign  of  dielectric 
antennas  of  these  shapes,  the  experimental  data 
must  necessarily  be  used.   It  will  latar  be  evident 
that  the  performance  of  rectangular  rods  is  very 
de 
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definitely  different  from  that  of  circular  rods. 
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CHAPTER  5. 

DESCRIPTION  OF  SOUIPMENT,  PROCEDURE,  AND  ET^ERIM^ITAL  DATA 

The  objective  of  the  experimental  v7ork  was  to 
determine  the  gain  on  axis  for  dielectric  antennas 
excited  In  the  dominant  mode,  using  optimum  length 
rods  in  each  case.   These  '.'/ere  measured  for  cylindrical 
rods  of  varying  cross  section  diameters,  and  for 
rectangular  rods  v;lth  varying  sizes  In  the  E  dimension, 
holding  the  H  dimension  constant  at  the  v;ldth  of 
the  standard  3  centimeter  vave  guide  (0.9  Inches). 

A  wavelength  (  Ao)  of  3.2  centimeters,  equivalent 

g 
to  a  freouency  of  9.37  x  10-^  cycles  per  second,  v;as 

used  for  the  reason  that  the  components  i^rere  of 

convenient  size  and  the  equipment  for  use  at  this 

frequency  vas  available.   The  equipment  iiaoBho^n  in 

Figures  10-14.   T>;o  dielectric  materials  v;ere  used 

for  each  antenna:  Polystyrene  (  £  -  ^.Sd   )  and  Teflon 

(^-J./O   ). 

A  standard  Klystron  oscillator  v:as  used  for  a 

transmitter,  but  square  wave  modulation  vras  Imposed 

on  the  carrier  frequency  to  facilitate  detection  in 

the  receiver.  The  transmitter  operated  from  a 
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regulated  power  supply  which  was  fed  from  the  60 
cycle  115  volt  pov;er  lines  through  a  voltage 
regulating  transformer.  The  receiver  consisted  of 
a  crystal  detector,  amplifier,  and  voltmeter.  As 
a  check  on  receiver  operation,  all  measurements 
were  repeated  usin'^;  a  Pickard  and  Burns  Bolometei^ 
Amplifier,  model  100,  Results  v:ere  in  substantial 
agreement  with  those  reported  and  \illl   not  be 
repeated. 

The  various  sizes  of  dielectrics  v;ere  fitted 
into  sections  of  v;ave  guide,  made  especially  for 
the  purpose,  tapered  dovm  to  standard  size  and 
fitted  V/'ith  a  standard  choke.   See  photographs.   The 
part  of  the  dielectric  which  fits  into  the  guide  was 
tapered  in  each  ca.se,  in  order  to  minimize  reflections 
and  reduce  standing  waves.   The  elimination  or 
correction  for  the  standing  waves  proved  to  be  the 
most  difficult  problem  to  overcome.   The  taper 
reduced  the  voltage  standing  wave  ratio  (VSv.Tl)  in 
most  cases  to  belov/  2.0.   In  every  case  this  v;as 
carefully  measured  and  a  correction  applied  to  obtain 
the  actual  transmitted  power,  using  the  v:ell  known  formula 

(193) 
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where  P+  Is  the  trannmltted  -oover 
Tm    Is   the  incident  pov;er 
r  Is  the  voltage  standin,^  wove  ratio  (V5"/R) 

Several  methods  were  used  to  make  the  measurements. 
To  be^in  v^lth,  the  dielectric  was  used  as  the 
transmitting^  antenna,  since  receiver  "oower  v;as 
inadecuate  to  accurately  measure  the  VS'iTR.   The  VSVJR 
was  measured  at  the  transmitter  using  a  standard 
standing  wave  Indicator.   Later  it  was  found  more 
precise  to  use  a  calibrated  attenuator  together 
with  a  slotted  section,  crj'^stal  detector,  am-JDllfier, 
and  voltmeter,  attenuating  as  necessary  to  equate 
the  voltage  from  the  high  voltage  point  to  the  lov/ 
voltage  point.  The  method  by  which  the  final 
measurements  were  made  used  this  equipment,  employing 
an  attenuator  which  had  been  accurately  calibrated 
by  the  National  Bureau  of  Standards. 

As  a  further  check  on  the  accuracy  of  the 
measured  gains,  the  various  dielectric  antennas 
v;ere  each  put  at  the  receiver  \/ith  a  special  matching 
section,  and  the  measurements  rer^eated.   For  each 
reading,  the  matching  section  v;as  adjusted  to 
maximize  readings  by  the  reduction  of  the  standing  vjaves, 
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Results  for  the  dielectric  antennas  up  to  one 
v.'avelength  in  diameter  or  in  the  S  dimension  vrere 
in  a?rreeraent  v/ith  those  measured  by  the  mrevious 
tv;o  methods.   G-ains  measured  for  larger  antennas 
were  somevrhat  less,  which  can  be  accounted  for  by 
the  fact  that  standing  waves  still  existed  to  an 
appreciable  amount  despite  the  matching  section. 

Some  slight  difficulty  v;as  encountered  vrith 
reflections  within  the  room.   These  \'ere  reduced 
as  much  as  r)os?ible  by  the  strategic  use  of  an 
3-bsorbing  screen.   As  a  check  on  errors  which  might 
be  introduced  due  to  reflections,  the  dists-nce  from 
the  transmitter  to  receiver,  normally  of  the  order 
of  one  hundred  vravelengths,  was  changed  in  ouarter 
wavelength  Increments,  noting  povrer  detected  at 
each  range.   If  a^-oreciable  reflections  v;ere 
interfering,  there  should  be  alternate  addition 
and  cancellation  of  -novrer  observed.   It  was  found 
that  detected  rjov/er  was  not  ar)preclably  different 
for  these  loositlons  a  quarter  v/avelength  apart, 
and  it  was  concluded  that  any  errors  from  reflections 
v;ere  of  little  significance. 

Before  mes-suring  the  gains  it  was  necessary 
first  to  find  the  optimum  length  of  the  dielectric. 
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Th.ls  v/as  done  experimentally  by  obGervln,^  relative 
power  transmitted,  first  vitli  the  dielectric  almost 
flush  v/ith  the  raouth  of  the  r/^ide,  and  then  out  l 
inch  and  so  on  with  lon.-jer  lenp;ths  using  increments 
of  about  -5  ^o   ,  Near  the  peak  of  the  transmission 
smaller  increments  v;ere  used  to  obtain  the  exact 
optimum  length.  Having  the  relative  nov/er  trans- 
mitted at  optimum  length,  the  test  antenna  was  then 
immediately  replaced  by  a  precalibrated  standard 
horn  and  the  T)0wers  compared.   This  enabled 
calculation  of  the  relative  gain  of  the  test  antenna 
with  respect  to  the  standard  horn. 

The  gain  of  the  standard  was  carefully 
measured  in  order  to  be  able  to  determine  the 
absolute  '^ain;  that  is,  the  ratio  of  the  peak  intensity 
of  the  given  antenna  to  the  intensity  of  an  isotropic 
radiator  when  both  are  radiating  the  same  total 
power.   In  terms  of  the  effective  area,  Ae  »  ^tie 


?'"  23 
gain,  G-,  has  been  shov.Ti "  '  -^  to  be 


47r/4 


e 


(194)     G    =   — -^ 

The  effective  area  Is  the  ratio  of  the  power 
received  at  the  antenna  to  the  -power  ner  unit  area 
in  the  incident  v/ave. 
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The  absolute  gain  of  the  standard  horn  v;as 
measured  by  several  of  the  methods  vrhich  have  been 
developed  and  Duolished  *   ' "'  .   One  of  these   will 
be  briefly  described  here,  -/hich  makes  use  of  tv;o 
identical  antennas.   The  standard  horn  i/as  attached 
to  the  transmitter  and  an  Identical  one  used  on  the 
receiver.   Under  these  conditions  it  is  easily  shov.Ti^ 
that 


(195)       p^   -   Ifr^' 


where  f^  is  the  received  Doi/er,  /^y-  is  the  transmitted 
TDOwer,  and  d  is  the  distance  between  the  tranem-itting 
and  receivini3  antennas.   To  determine  the  ratio  of 
received  "oov/er  to  total  transmitted  pov;er,  the 
transmitter  and  receiver  were  coupled  directly 
together,  exc-f>pt  for  a  calibrated  attenuator  betv/een. 
v.- 1th  a  large  but  accurately  measured  amount  of 
attenuation  inserted,  a  particular  receiver  voltage 
was  noted.   Then  the  two  antennas  v.'ere  attached 
respectively  to  the  trp-nsmitter  and  receiver  and 
separated  by  an  accurately  measured  distance 
(several  feet),  and  the  attenuation  of  the  transmitter 
reduced  until  the  same  receiver  voltage  v/as  obtained. 
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The  difference  in  the  amount  of  attenuation  In  the 
two  cases  gave  the  ratio  of  Pp  to  ?„  In  decibels, 
from  v;hlch  G  could  easily  be  calculated. 


Presentation  of  Data 

Three  alternatives  in  the  presentation  of 
experimental  data  have  been  considered,  and  these 
v;ill  be  briefly  mentioned  to  show  the  reasons  for 
the  selection  of  the  one  /^iven. 

(a)  Report  the  relative  r-;ain  of  the  dielectric 
antenna  compared  to  the  ^ain  of  the  feeding  v/ave 
guide  v;ith  the  dielectric  removed.   This  v/ould 
permit  a  direct  comparison  of  the  pov;er  on  axis 
with  and  v:lthout  the  dielectric  vrhich  v;ould  ohov; 
clearly  the  effect  of  the  dielectric  in  directing 
the  beam.   However  this  method  was  excluded,  because 
the  experimental  data  v;ould  not  be  reproducible 
unless  the  feed  structures  were  identical.   G-aln 

of  open  mouth  feeds  are  not  precisely  rredic table 
and  there  would  be  no  way  to  obtain  absolute 
gain  oi"  effective  area  from  these  data, 

(b)  Report  the  absolute  gain  of  the  dielectric 
antenna;  that  in  the  oower  on  axis  relative  to  an 
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isotropic  oscillator  emittln,'^  the  same  total  pov/er. 
This  is  better,  but  here  a^ain  variations  in  the 
metallic  feed  structures  v;o.uld  cause  variation  in 
the  absolute  gain.   Furthermore,  as  shown  by 
equation  (194)  above,  the  scales  v;ould  have  an 
artificial  multiplier  of  4T. 

•(c)   Report  the  effective  cross  section  area 
relative  to  the  "ohysical  cross  section  area.   This 
is  the  method  chosen  because:  (l)  this  presentation 
eliminates  the  undesired  effect  of  variation  in 
feed  structures;  (2)  as  cross  section  areas  are 
allowed  to  grov/,  the  effective  cross  section  areas 
also  »^rov,  and  these  tend  to  become  eaual.   In  other 
v.'ords,  the  captured  energy  tends  tov/ard  that 
physically  interce-nted  by  tliis  area.   Therefore, 
a  plot  of  effective  area  versus  physical  area  v/ill 
verge  tcard  a  straight  line  irlth   unit  slope  passing 
throu.jh  the  origin  regardless  of  the  dielectric 
used;  and  (3)  an  extremely  useful  piece  of  information 
concernirvg  any  antenna  is  its  effective  cross 
section  area.   From  this  figure  it  is  a  simple  and 
direct  calculation  to  obtain  the  total  energy 
gathered  by  the  antenna  in  a  field  of  a  given  strength. 

This  third  method  of  presentation,  in  addition 
to  providing  the  engineer  the  very  useful  information 
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of  maximum  capture  area  for  various  dielectrics 
of  various  sizes,'  illustrates  the  variation  of 
capture  area  from  the  actual  area  in  the  sizes 
of  the  order  of  a  v/avelength,  and  shov/s  hov;  these 
should  converge  to  the  45  line  as  dimensions 
increase  tov/ard  o-otical  sizes. 


Choice  of  Scales 

Consideration  was  given  to  reporting  dimensions 
directly  in  centimeters  or  inches.   This  v;ould 
enable  an  en/^ineer,  employing  v/avelengths  the  same 
as  used  here  (3.2  cm)  to  use  the  data  directly. 
However,  this  was  rejected  because  the  value  of 
these  data  is  not  limited  to  any  particular  \7avelength. 
Dimensions  are  therefore  given  in  non  dimensional 
form;  that  is,  lengths  in  wavelengths  and  cross 
section  areas  in  square  wavelengthc .   Although  this 
requires  conversion  of  practical  units  to  this  system 
to  enable  usage  of  the  data,  it  has  the  advantage 
that  the  data  are  applicable  to  a  very  wide  band 
of  wavelengths. 

A  plot  of  effective  area  versus  actual  cross 
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section  area  for  the  cylindrical  dielectric  antennas 
measured  is  ^iven  in  Picture  15.   Flp^ure  16  is  the 
same  plot  for  rectangular  cross  section  antennas. 
Figures  17  and  18  are  also  enclosed  for  completeness 
(giving  the  absolute  gains  of  cylinders  and 
rectangular  cross  sections,  respectively,  versus 
cross  section  area. 

It  iiTill  be  noted  that  measurements  v;ere  not 
obtained  for  cylinders  of  lar^e  enough  diameter 
to  shovr  convergence  of  the  effective  area  tov.'ard 
the  actual  area.   In  the  larger  diameter  circular 
cylinders,  the  mode  tended  to  shift  to  the  /?}  =  o 
mode  v:hich  gives  a  null  on  axis.   No  v.-ay  to  nrevent 
this  v^as  devised.   There  iras  also  a  practical  linit 
to  the  measurements  on  the  low  side,  ^-rith  the  time 
and  facilities  a,vailable.   For  small  cross  section 
rods  thd  optimum  length  becomes  unmanageably  lonp^, 
and  the  thin  rods  tend  to  sag.   Furthermore,  these 
thin  long  rods  require  a  long  range  in  order  to 
put  the  field  meter  in  the  far  zone. 

The  present  technique  does  not  permit  extension 
of  the  range  of  measurement  for  circular  rods, 
and  the  possibility  of  existence  of  TM  modes  in 
these  dielectrics  is  a  disturbing  factor.   The 
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apparent  irregularity  of  the  polystyrene  rods  may 
be  caused  by  this  r)henomenon.  In  order  to  reduce 
the  probability  of  spurious  modes  in  the  excitin,^ 
radiator  for  the  larger  areas,  conical  horns  v/ere 
also  tested  v^ith  polystyrene  cylinders,  v/lth  the 
following  results 


Horn  diameter 


Ati 


Horn 
Physical  Area 

Antenna 

Capture  Area 

2.5  ^-" 

3.T8/U 

5.25  A^ 

5.24  ;ir- 

2t 


3i" 


Here  the  trend  of  antenna  effective  area  tovrard 
actual  area  is  noted.   The  horns  v/ere  far  from 
optimum  from  the  point  of  view  of  radiation  pattern, 
however,  and  the  lens  antennas  were  alno  probably 
not  optimum.   Therefore,  convergence  of  optimum 
cylinders  tovrard  the  expected  asymptote  v;ould  not 
be  closely  realized  until  the  experiment  reached 
much  larger  areas  than  the  above  figures  indicate. 

In  the  rectmgular  rods,  no  tendency  of  the 
mode  to  shift  vjas  observed,  but  other  practical 
difficulties  prevented  extension  of  the  rapge  of 
sizes.   In  order  to  prevent  cut  off  in  the  antenna 
feed,  v;hen  using  small  sizes,  the  H  dimension  was 
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set  at  a  constant  value  of  .71  >^o  .  Small  areas, 

then,  v/ould  have  required  very  lon,^  and  thin  sections  J^N 

which  would  have  caused  the  same  difficulties  mentioned 

above  for  circular  rods.   The  larger  cross  section 

areas,  v;here  the  E   dimension  v^as  large  compared 

to  the  constant  H  dimension,  gave  radiation  patterns 

which  were  relatively  narrow  in  the  -nlane  -narallel 

to  the  E  dimension,  and  relatively  broad  in  the 

plane  -oarallel  to  the  H  dimension,  thus  giving 

reduced  gain  on  axis.   It  v;as  found  that  if  the 

cross  section  area  had  heen   made  more  nearly  snuare, 

a  radiation  pattern  more  nearly  circular  vjould  have 

been  produced,  and  th3reby  increasing  gain  on 

axis.   The  gain  on  axis,  being  a  function  of  the 

radiation  pattern,  then,  is  der>endent  on  the  shape 

of  the  cross  section.   For  small  cross  sections  this 

effect  is  small,  but  for  cross  sections  of  a  vjave- 

length  or  larger  the  effect  is  very  Tnronounced. 

A  v:ay  to  avoid  this  vjould  be  to  properly  shape 

the  front  surface  of  the  dielectric,  or  to  to/oer 

the  rod.   Neither  of  these  were  attempted  here. 

The  end  surfaces  of  all  the  antennas  tested  v/ere 

plane. 

It  v-'as  evident  that  the  larger  rectangular 
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dielectrics  v;ould  /^ive  higher  Dov^ers  on  axis  if  the 
cross  section  \vas  apDroximately  square.   Therefore, 
as  the  v/avelength  region  is  approached  in  Figure  16, 
the  larger  physical  cross  section  antennas  appear 
at  a  disadvantage.   Perlmps  a  better  comparison 
in  this  region  vrould  have  been  to  compare  the  antenna 
capture  area  to  the  capture  area  measured  for  the 
same  metallic  feed  v;ith  the  dielectric  flush  v;ith 
the  metallic  aperture  (zero  length  antenna).   This 
method  of  reporting  would  have  brought  the  larger 
areas  more  in  accord  vith  the  following  table, 
measured  only  for  nolystyrene. 


Feed  Korn    Physical  Area  Ganture 


Car)tnre  Area 


Dimensions 

in   Ao^ 

Area 
Feed  Horn 

of 

Antenna 

.9  X  1.8" 

1.02 

0.64 

1.03 

1^   X   1^" 

1.42 

1.0 

1.7 

2%-  X  2^4" 

4.76 

4.75 

5.3 

This  table  shov;s  the  trend,  predicted 
previously,  of  the  capture  area  of  the  optimum 
length  antenna  to  approach  th;^  physical  cross 
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section  as  the  latter  beoomes  many  v;avelen^ths  in 
each  direction.   The  shaping  of  these  larger 
antennas  v;as  only  attempted  approximately  and 
somewhat  larger  capture  areas  v:ould  be  expected 
if  optimum  shapes  had  been  employed. 
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CHAPTER  4. 

DISCUSSION  OF  RESULTS  AND  ATTSI^TPTED  CORRSLATION  V;iTH 

EXISTING  THEORIES 

The  curves  of  effective  area  versus  cross  section 
area.  Figures  15  and  16  of  Chapter  3»  show  that 
considerable  increases  in  effective  areas  over 
actual  areas  are  obtainable  if  the  physical  areas 
are  small  (less  than  a  square  v:avelength) .   These 
effective  areas  are  auite  high  for  small  physical 
areas  and  reduce  toward  the  physical  areas  as  sizes 
become  lo,rger.  The  trend  is  apparent  for  the 
rectangular  dielectrics,  but  could  not  be  demonstrated 
with  circular  cross  section  because  of  shifting 
of  the  Siode  v.rithin  the  dielectric  despite  excitation 
by  the  ii  =  1  mode. 

As  previously  pointed  out,  the  design  of  the 
dielectric  antenna  above  a  i;avelength  in  cross 
section  area  is  still  not  understood.   The  gain  on 
axis,  on  v/hich  the  experimental  results  are  based, 
is  a  function  of  the  sha-oe  of  the  total  -oattern. 
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The  gain  is  measured  by  samplinr;  this  pattern  on 
the  prolonged  axis  of  the  dielectric.   This  STDread 
in  iiattem  is  -oarticularly  apparent  in  these  data 
in  the  sizes  Immediately  above  a  souare  wavelength. 
In  this  region  ray  optics  play  an  inrortant  role; 
and  yet  customary  treatmonts  by  ray  optics  fail 
because  of  the  considerable  phase  gradients  due 
to  the  adjacent  metallic  structures  supporting  the 
radiator.   In  fact,  such  a  treatment  v;ould  involve 
serious  difficulties.   In  the  present  project,  the 
problem  would  be  further  complicated  by  the  absence 
of  an  accurate  means  to  det'^rmine  the  complete  vrave 
characteristics  of  each  of  the  various  feed  horns 
used.   Proper  lens  shaping  v:ould  necessarily  have 
to  be  based  on  a  knowledge  of  the  ^hase  front  of 
the  horn  to  be  lens^d.   However,  to  shov:  the  trend 
beyond  the  one  ■vravelength  region  tov/ard  the  expected 
45  asymptote  a  few  measurements  were  attemDted 
in  V7hich  the  face  of  the  dielectric  'was  shaped  by 
means  of  trying  many  combinations  of  various  si-.e 
blocks  so  as  to  obtain  a  stepped  profile.   The 
best  gains  v;ere  found  in  these  larger  shapes  when 
the  Drofiles  resembled  o-otical  lenses  in  physical 
shane,  and  the  larger  the  area  the  more  nearly  the 
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true  optical  shape  is  reouired  to  optirdze  gain. 
This  is  in  contrast  to  smaller  cross  section  areas, 
where  sha-nin^^  of  the  front-  face  had  little  effect, 
and  exnerirasntal  measurements  v^ere  very  consistent. 
Evidently  the  phase  gradients  vrithin  the  smaller 
dielectrics  play  a  much  less  significant  role. 

Because  of  the  fact  that  these  data  are  based 
on  unoptimized  profiles,  it  is  apparent  that  gains  of 
optimized  sha-oes  would  "be  liigher  than  here  reoorted, 
and  convergence  to  the  45  asymptote  vrould  occur  much 
further  out.   It  is  exoected  that  when  better  theories 
are  available  for  the  design  of  lenses  in  the  v/avelength 
region,  curves  similar  to  those  of  Figure  16  could  be 
expected  to  converge  monotonlcally  to  a  line  of  unit 
slope,  and  to  reach  the  asymptote  only  after  dimensions 
are  such  that  the  areas  are  ms.ny  square  wavelengths. 

One  experiment  was  attempted  some  years  ago  to 
"lens"  a  30  inch  parabolic  aperture  on  X  band,  and  it 
was  not  found  possible  to  detect  any  increase  in 
effective  area.   As  the  aperture  grovrs,  optical  lenses 
are  obtained,  where  effective  area  is  strictly  proportional 
to  actual  area.   In  a  recent  report'^  of  experiments 
sponsored  by  the  U.S.  Array  Signal  Corps,  tests  v:ere 
revealed  on  a  spherical  Lucite  lens  of  cross 
section  area  985  square  wavelengths.   G-ains 
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reported,  after  correction  for  the  losses  occurring 
in  Lucite,  indicate  an  effective  area  of  635  square 
wavelength-s.   The  apparent  discrepancy  is  due  to 
the  fairly  narrow  illumination  used  in  this 
experiment.   The  illumination  duplicated  the  usual 
radar  dish  illuminations  and  resulted  in  the 
expected  aperture  utilization  factor  of  65^ 
commonly  employed  in  radar  dish  antennas.   This  shows 
that  effective  areas  actually  do  converge  to  the 
physical  areas  utilized  in  large  sizes. 

Converging  lenses  have  been  made  for  the 
vravelengths  considered  herein.  V/hile  these  appear 
to  give  a  tremendous  pov/er  at  their  designed  focal 
point,  a  more  thorough  discussion  v/ould  show  that 
they  merely  operate  as  converging  devices,  obeying 
to  a  first  approximation  the  half  povrer  law 


4-  =  ^   TDOwer  angle  of  radiation 


They  do  not  increase  the  effective  area  of  the 
aperture. 


Mars ten* 3  Theory 

A  very  excellent  theoretical  treatment  of  the 
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dielectric  antenna  of  circular  cross  section  is 
given  by  Marsten  in  a  recent  report   published 
by  the  U.S.  Naval  Research  Laboratory.   A  circular 
dielectric  rod  is  excited  in  any  sinr^le  mode,  and 
the  equations  used  for  the  v.-aves  are  equivalent  to 
those  given  in  Chapters  1  and  2.   From  these  are 
calculated  equivalent  electric  and  magnetic  current 
distributions  over  the  surfaces,  in  accordance  v/ith 
Schelkunoff 's  Equivalence  Theorems   .   Then  the 
radiation  field  is  computed  usin^^  the  method  of 
Schelkunoff,  first  for  the  field  arising  from  the 
equivalent  sources  o,t  the  sices  of  the  rod,  £^' 
and  then  for  the  equivalent  sources  at  the  end  of 
the  rod,  £   ,   Gerta.in  factors  common  to  both 
expressions  are  dropTDed,  since  only  the  relative 
field  is    ^  of  interest,  and  the  final  expressions 
for  the  9  components  are  given  as  follov/s 


(196) 


/l^§cose-uJ^^,(L) 
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(196)  continued 


!i±L-sL^iLii  jjoj^fr)  -  ^Jj^)j:a)] 


The  expressions  for  the  other  com^^onents  of  the 
field  will  not  be  repeated,  since  only  the  pov:er 
on  axis  is  of  interest  here. 

The  notation  is  the  same  as  given  herein, 
except  as  follov;s: 

Jl   _  —  _  len/^th  of  rod.   Since  optir.un  length  fi'&j 


is  here  considered,  substitute  ^° 


X  -    t(x  c^sa  -  -^) 


TT 


On  axis  Cos 6"  =  1  and  for  optimum  length  X-  -  ^ 

L  ^  j0^  ^//7  6»  On  axis  i^  0 

A  and  B  are  functions  which  are  always  finite 
and  non-zero 
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S 


which  UT^on  sirarllflcction  a.nd  substitution  from 
equations  (I87)  and  (194)  becomes  in  our  notation 

By  actually  substituting  numbers,  It  is  found 
that  6   has  a  plus  value  for  all  dimensions  under 
consideration,  and  so  the  ;'^ositlve  sir^,n  should  be 
used  after  extractin.^  square  root, 

Assuminr;;  n  =1,  e  =  o  .  i  ^  -3^  ,  ^^^  p  ^0 
the  above  eauatlons  reduce  to 


(197) 


(198) 


e;  =  -^iiff  j^vp)[n^.^o  iE^^jA. 
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A  plot  of  effective  area  for  optimum  length  comouted 

as  above  Is  given  in  Figure  20,  vhere  the  abscissa 

is  taken  as  the  physical  area. 

The  gains  calculated  by  Marsten's  theory  show 

the  same  general  trend  as  the  experimental  data 

of  Chapter  3  show.   The  principal  exceptions  are, 

first,  that  all  are  much  higher  in  value  than  could 

be  found  exT.)erimen tally,  and,  second,  it  is  noted 

that  the  theory  does  not  hold  at  all  when  dimensions 

are  small  enough  so  that  n^^  approaches  1.   Equation 

(202),  the  expression  for  the  contribution  from  the 

sides,  has  a  denominator,  ]rf^<*.~0     which  approaches 

zero  as  n   aDDroaches  1  v:ith  decreasing  diameters. 
a   ' ' 

The  term  -J^ffj  in  the  numerator  also  approaches  zero 
as  the  diameter  goes  to  zero,  but  this  occurs  much 
more  slowly.   The  result  is  that  pov/ers  on  axis  as 
indicated  become  extremely  high  for  small  rods.   These 
enormous  pov;ers  do  not  conform  to  experimental 
observations. 

The  question  of  identical  phases  for  the 
radiation  from  the  side  walls  and  from  the  ends 
does  not  appear  to  be  correct  from  a  physical  viewpoijit. 
The  -ohase  v;-ithln  an  optimum  length  rod  emerges  one 
half  period  later  than  the  v/ave  travelling  the  same 
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distance  in  free  s-oace.  Alon^r^  the  rod  the  phase 
relative  to  the  end  must  certainly  lead  by  something 
less  than  a  half  period  but  raore  than  zero.   If, 
In  the  theory,  some  adjustment  to  these  r)hases 
were  made,  it  Is  possible  that  the  theory  might 
fit  the  experimental  flndin^js  more  closely. 

Another  interesting  point  is  noted  from 
equations  (196)  and  (197).   As  the  diameter  becomes 
very  large  p  approaches  zero.  (  p-^  a \/\^f--^ 

and  at  large  diameters  h  approaches  k-j^) .   If  n  «=1, 
then  J^  (p)  =  Oj  9-nd  the  right  sides  of  both  equations 
vanish.  V/e  must  conclude  from  this,  If /I  theory 
holds  for  large  sizes,  that  no  radiation  could  occur 
from  an  antenna  of  very  large  diameter  if  the 
optimum  length  is  determined  as  •calcula1;.ed  above. 


vqikes*  Theory 

An  explanation  for  the  peculla,r  behavior  of 
gains  in  the  region  of  areas  of  a  fraction  of  a 

square  v.^avelength  is  given  In  a  theory  proposed  by 

24 

V/llkes   on  the  Interaction  of  Radiant  Energy.   It 

has  been  shown  experimentally  that  the  effective 


104. 


areas  in  this  region  are  considerably  greater  than 
the  -Dhysical  cross  section  areas,  although  pov:ers 
are  definitely  limited.   Presently  published  theories 
on  dielectric  antennae  of  circular  cross  section 
v/ould  allov/  an  unlimited  growth  of  capture  area 
v;ith  diminishing  cross  section  areas. 

l/ilkes,  however,  has  shoim  experimentally  that 
the  -oower  grov;th  is  in  reality  not  unlimited,  but 
that  the  captui^e  area  for  small  cross  sections  Is 
limited.  A  simplified  qualitative  outline  of  his 
approach  is  given  below  to  indicate  the  method  by 
\-/hich  he  arrives  at  a  solution  of  this  hitherto 
unaccountable  phenomenon. 

As  the  propagation  dovm  a  very  thin  dielectric 
rod  is  a  plane  wave,  the  phase  and  group  or  photon 
velocities  are  assumed  by  liim  to  be  identical.  By 
an  ingenious  technioue,  he  has  been  able  to  measure 
the  phase  velocities  in  extremely  short  lengths  of 
dielectric  and  to  compute  the  velocities  down  to 
a  zero  length.   In  these  thin  dielectrics  the 
energy  starts  at  its  critical  velocity,  -^     for 
very  short  lengths  and  aporoaches  c  for  increasing 
lengths.   Knowing  this  rate  of  change  of  velocity 
it  is  possible  to  conceive  of  the  photons  following 
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describable  ballistic  trajectories.   This,  hovrever, 
requires  the  concept  of  an  applied  force  to  photons 
of  finite  mass  v/hich  exhibit  properties  of  inertia. 
Since  the  photon  mass  is  at  present  unl^no^.n,  except 
for  de  Bro-^lie's  controversial  findings"  ,  direct 
use  of  ballistic  principles  to  calculate  this 
force  is  not  possible.   Nevertheless,  with  this 
concept  in  mind,  a  consideration  of  the  nature 
of  the  forces  brought  to  bear  on  the  photon  reveal 
a  behavior  closely  compatible  -rith  exfjerimental 
observations. 

Postulatin>s  the  existence  of  forces  betv;een 
the  internal  and  external  fields,  v;hich  v;e  have 
Just  noted,  it  follovrs  that  there  is  a  transfer 
of  work  or  energy  from  one  of  these  fields  to  the 
other.   The  medium  of  energy  exchange  is  photons 
of  a  frequency  determined  by  the  excitation, 

Wilkes  has  found  experimental  evidence  of 
the  existence  of  these  forces;  in  fact  his 
calculations  and  experiments  support  the  existence 
of  the  force  as  composed  of  tv/o  parts.   One  part 
Is  a  pressure  action  proportional  to  the  difference 
in  velocities  of  the  tv:o  'fields.   The  other  is  a 
function  of  the  total  ohase  lag  and  a-nr^ears 
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gravitational  in  nature.   The  ma5nitude  of  this 
action  is  foumd  eoual  to  the  radiant  pressure  and 
parallels  in  this  respect  Mandelker's  theory  '   of 
gravitation. 

From  Quantum  Electrodynamics,  vrilkes  finds 
that  the  quanta  in  the  accelerated  dielectric 
fields  become  negative  and  -nresent  holes.   He  then 
allows  these  holes  to  be  filled  in  to  the  limit 
allov/ed  by  the  statistical  lav/s  of  a  degenerate 
gas,  and  by  means  of  these  calculations  arrives 
at  the  tv/o  saturation  joints  shown  on  Figure  21  as 
the  loi/est  cross  sections  reported.   These  points 
have  been  confirmed  by  his  experiments.   The  other 
points  of  Figure  21,  namely  those  for  cross  sections 
greater  than  0,2  square  wavelengths,  are  8. 
repetition  of  the  data  obtained  in  Chapter  3. 
Obviously,  if  the  physical  area  drops  to  zero,  the 
capture  area  vrill  also  go  to  zero,  and  this 
seeminp;ly  trivial  statement  permits  the  completion 
of  the  curves  of  Figure  21. 
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Conclusion 

It  is  believed  that  no  theory  based  purely  on 
Maxvrell's  eouations  vrill  suffice  to  explain  the 
neculiar  behavior  of  gains  In  the  very  small  cross 
section  region.  As  stated  early  in  this  oaper, 
Max^fell's  eouations  hold  vrith  certainty  only  for 
large  scale  electromagnetic  phenomena,  and  r:hen 
dimensions  are  reduced  to  fractions  of  a  wavelength, 
they  cannot  be  relied  upon.   The  true  explanation  of 
very  small  dielectric  antennas  must  come  from  quantum 
electrodynamics,  as  proposed  by  V/ilkes,  and  the 
data  obtained  in  Chapter  3  appear  to  be  in  close 
agreement  with  this  theory. 

As  dimensions  grow  larger,  effective  areas 
of  dielectric  antennas,  if  shaped  for  optimum  gain, 
remain  greater  than  the  physical  cross  sections,  but 
will  tend  monotonically  toward  these  areas.   After 
areas  ha.ve  reached  several  square  wavelen^gths, 
effective  area  of  optimum  shaped  dielectrics  can  be 
closely  approximated  by  the  actual  physical  area. 
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GLOSSARY  OF  3Y130LS 

E  31ectric  Field  Intensity  vector 

H  Magnetic  Field  Intensity  vector 

D  Slectric  Displacement  vector 

B  Ma^^netic  Induction  vector 

J  Current  Density  vector 

f         Charge  density  at  any  noint.   Also  used  as  a 
general  argument  of  a  Bessel  function 

P    Electric  -oolarization  vector 

M    Magnetic  polarization  vector 

g^    Electric  Inductive  capacity  of  free  space 
defined  by    D=  ^E 

yC{^        Magnetic  Inductive  capacity  of  free  s^^ace 
defined  "oy  6-MoH 

X^   Magnetic  susceptibility     M^  X^// 

"Xg   Electric  susceptibility     P  =  T^e^ 

£    Relative  dielectric  constant  -  -f- 

Co 


K^  Permeability  -  .. 

A  Vector  Potential  of  the  Field 

(|)  Scalar  Potential  of  the  Field 

II  Hertz  vectorj-Slectric  Polarization  Potential 

II  Hertz  vector.  Magnetic  Polarization  Potential 

^  Conductivity 
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W  Surface  density  of  charge^ or  charge  per  unit  area 

I  Conduction  current 

1  Square  root  of  minus  one 

h  Longitudinal  propagation  factor  «  ^ 

t  Time 

e  Base  of  natural  logarithms 

-J^(f]  Bessel  Function  of  order  n 

f  An  argument  of  a  Bessel  function 

7  (p)   Circular  cylinder  function,  or  Bessel  function, 
^^^■*^   of  order  n 

W^[^]   Bessel  function  of  the  second  kind 

H^(f}Hankel   function  of   the  first  kind 

L(?)  Modified  Bessel  function 

k    Intrinsic  ^hase  constant  -  yMiuy"  ^'X^e^ 

0)         Angular  velocity  =  freouency/^r 

X         A  propagation  constant  for  dielectric  in  the 
case  of  a  circular  cylinder  in  air  =  ]/j^~-'^ 

y,^      A  propagation  constant  for  air  in  the  case  of 
a  circular  cylinder  in  air    l/^J^^ 

p    A  propagation  constant  for  dielectric  =  x,<j? 

V    A  propagation  constant  for  air  in  the  case  of 
a  circular  cylinder  in  air  =  ?(i  i^ 

a    Radius  of  a  circular  cylinder  of    dielectric 

q    A  propagation  constant  for  air  =  -Uyir  -   -.<Xj«.  = -xa.\r^~J^ 
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n    Apparent  index  of  refraction  for  dielectric  rod 


f    Abbreviation  for 


I 


g         Abbreviation  for        ^,y  ^f\ 

Ac        Free   SDace  vjave length 

7l,         V.^avelength  of   propagation  in  a  dielectric  a^ntenna 

c  Velocity  of  Dropar^ation  in  free   space 

J^o        Optimiim  length  for  dielectric   to  maximize 
on-axis  radiation 
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